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CHAPTER  I 
INTRODUCTION 


1.1  Background 

There  is  a  strong  motivation  for  the  study  of  the  vibrations 
of  cylindrical  shells.  They  are  the  prototypes  of  such  practical  struc¬ 
tures  as  pipelines,  containment  vessels,  undersea  vehicles,  aircraft, 
and  machine  housings.  Insight  into  the  dynamical  response  of  cylindri¬ 
cal  shells  to  various  exciting  forces  is  essential  to  the  understanding 
of  the  sound  radiation  from,  or  the  vibration  isolation  between  com¬ 
ponents  of  many  practical  structures. 

Compared  to  structural  elements  such  as  beams  and  plates,  the 
shell  vibrates  in  an  extremely  complex  manner.  Because  of  the  shell's 
curvature,  it  is  impossible  to  separate  the  extensional  or  membrane-like 
motion  from  the  bending  or  plate-like  motion.  With  coupling  between  ex¬ 
tension  and  bending,  the  development  of  an  equation  of  motion  is  a  for¬ 
midable  task.  Of  all  the  possible  shell  structures,  only  the  circular 
cylindrical  and  spherical  shell  are  mathematically  tractable.  .An  ex¬ 
cellent  and  self-contained  derivation  of  the  equations  governing  the  vi¬ 
bration  of  shells  is  presented  in  Reference  [1]. 

Much  of  the  early  literature  on  cylindrical  shells  was  di¬ 
rected  toward  the  development  and  solution  of  various  shell  equations  of 
motion  to  determine  natural  frequencies  and  mode  shapes  for  various 
boundary  conditions.  The  development  of  shell  theories  and  their 


1 


CHAPTER  I 
INTRODUCTION 


1.1  Background 

There  is  a  strong  motivation  for  the  study  of  the  vibrations 
of  cylindrical  shells.  They  are  the  prototypes  of  such  practical  struc¬ 
tures  as  pipelines,  containment  vessels,  undersea  vehicles,  aircraft, 
and  machine  housings.  Insight  into  the  dynamical  response  of  cylindri¬ 
cal  shells  to  various  exciting  forces  is  essential  to  the  understanding 
of  the  sound  radiation  from,  or  the  vibration  isolation  between  com¬ 
ponents  of  many  practical  structures. 

Compared  to  structural  elements  such  as  beams  and  plates,  the 
shell  vibrates  in  an  extremely  complex  manner.  Because  of  the  shell's 
curvature,  it  is  impossible  to  separate  the  extensional  or  membrane-like 
motion  from  the  bending  or  plate-like  motion.  With  coupling  between  ex¬ 
tension  and  bending,  the  development  of  an  equation  of  motion  is  a  for¬ 
midable  task.  Of  all  the  possible  shell  structures,  only  the  circular 
cylindrical  and  spherical  shell  are  mathematically  tractable.  .An  ex¬ 
cellent  and  self-contained  derivation  of  the  equations  governing  the  vi¬ 
bration  of  shells  is  presented  in  Reference  [1], 

Much  of  the  early  literature  on  cylindrical  shells  was  di¬ 
rected  toward  the  development  and  solution  of  various  shell  equations  of 
motion  to  determine  natural  frequencies  and  mode  shapes  for  various 
boundary  conditions.  The  development  of  shell  theories  and  their 


2 


associated  equations  of  motion  reached  a  high  degree  of  accuracy  (and 
complexity)  by  the  late  1950's.  Reference  [2]  provides  an  excellent 
overview  of  the  historical  development  of  various  thin  shell  theories. 

Of  particular  importance  in  the  development  of  modern  shell  equations  of 
motion  is  the  paper  by  Arnold  and  Warburton  [3] .  This  paper  quantified 
the  relationship  between  bending  and  extensional  energies  of  vibration 
for  the  cylindrical  shell. 

The  profusion  of  thin  shell  theories  led  Greenspon  [4]  to 
evaluate  the  range  of  applicability  of  these  theories  in  comparison  with 
exact  thick  shell  theories.  Greenspon  concluded  that  the  membrane  the¬ 
ory  of  shells  is  accurate  for  predicting  frequencies  and  displacement 
ratios  of  cylinders  with  appreciable  thickness. 

An  excellent  reference  for  the  researcher  and  engineer  is  the 
monograph  prepared  by  Leissa  [5].  This  monograph  summarises  over  500 
papers,  published  before  1977,  on  cylindrical  shells.  Chapter  I  provides 
a  very  detailed  discussion  of  the  differences  between  the  various  shell 
theories.  One  can  generally  conclude  that  the  differences  between  the 
shell  equations  arises  from  small  differences  in  the  strain-displacement 
relationships,  and  these  differences  have  little  significance  in  most 
practical  applications. 

The  development  of  the  shell  theories  has  produced  a  good  un¬ 
derstanding  of  the  behavior  of  vibrating  shells.  When  it  comes  to  prac¬ 
tical  applications,  however,  these  methods  do  not  go  very  far.  In  prac¬ 
tice,  we  are  usually  concerned  with  the  response  of  a  vibrating  system 
to  a  given  excitation  as  a  function  of  frequency.  Specifically,  we  are 
interested  in  the  resulting  velocity  on  the  surface  of  a  vibrator  for  a 
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given  force  distribution.  This  information  is  embodied  in  the  mechani¬ 
cal  impedance,  or  its  reciprocal,  the  admittance.  The  simplest  force 
distribution  is  a  point  force.  Because  the  response  due  to  an  arbitrary 
forcing  function  can  be  generalized  from  that  of  a  point  force,  we  are 
most  interested  in  the  response  to  point  excitation. 

Although  there  is  a  wealth  of  literature  available  concerning 
the  shell  equations,  resonance  frequencies,  and  mode  shapes,  there  are 
relatively  few  papers  which  address  the  mechanical  impedance  of  cylin¬ 
drical  shells. 

Franken  [6]  developed  an  expression  for  the  impedance  of  non- 
axisymmetric  modes  to  point  excitation  below  the  radial  resonance  fre¬ 
quency.  The  ring-stiffened  shell  problem  was  set  up  conceptually,  but 
was  not  solved  analytically.  The  models  developed  in  this  thesis  bear 
a  strong  resemblance  to  the  circuit  representation  in  Franken' s  paper. 

The  model  suggested  by  Franken  provides  a  good  estimate  of  the  ring- 
shell  response  at  low  frequencies,  but  at  high  frequencies  the  model 
must  be  modified. 

Based  on  approximate  shell  equations,  Heckl  [7]  developed  sim¬ 
ple  formulas  for  the  resonance  frequencies,  modal  densities,  and  driving- 
point  impedance  of  cylindrical  shells.  Although  his  results  are  based  on 
many  simplifications,  Heckl  clearly  demonstrates  that  the  simplifications 
are  more  than  accurate  enough  for  practical  vibration  control  work.  The 
simplified  frequency  equations  which  Heckl  presents  are  used  in  this  the¬ 
sis  in  deriving  the  density  of  resonances  for  a  cylindrical  shell. 
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Palladino  and  Neubert  [8]  demonstrated  the  use  of  mode-sum 
techniques  for  the  prediction  of  the  admittance  of  long  pipe-like  cylin¬ 
drical  shells  to  point  excitation.  Their  paper  also  provides  a  good 
discussion  of  the  anomalies  that  occur  in  transfer  admittance  data. 

Because  of  the  complex  nature  of  vibrating  structures,  it  is 
desirable  to  develop  simplified  methods  for  the  estimation  of  gross  re¬ 
sponse.  In  practice,  we  are  usually  concerned  with  only  the  mean  re¬ 
sponse  over  some  frequency  band.  Quite  often  the  excitation  can  only  be 
described  in  terms  of  broad  frequency  spectra,  and  in  general,  one  would 
be  interested  in  only  the  broad  frequency  response. 

Skudrzyk  [9,  10,  11]  has  developed  expressions  for  the  geo¬ 
metric  mean  of  the  admittance  of  the  most  important  vibrators:  rods, 
beams,  rings,  plates,  shells,  and  membranes.  Skudrzyk’s  Mean-Value  Ad¬ 
mittance  Method  predicts  the  mean  line  through  the  logarithmically  re¬ 
corded  frequency  response  curve  of  a  complex  vibrator,  the  height  of  the 
resonance  peaks,  and  the  depth  of  the  antiresonances.  The  results  ob¬ 
tained  using  this  method  provide  a  good  approximation  with  which  to  eval¬ 
uate  the  design  of  many  complex  structures  consisting  of  the  fundamental 
components  of  plates,  shells,  beams,  and  rings.  One  of  the  most  im¬ 
portant  features  of  this  method  is  that,  in  contrast  to  other  methods 
such  as  Statistical  Energy  Analysis  and  Finite  Element  Analysis,  the 
Mean-Value  Admittance  Method  is  applicable  at  low,  as  well  as  high  fre¬ 
quencies. 

1.2  Purpose 

The  purpose  of  this  thesis  is  to  demonstrate  a  practical 
method  of  predicting  the  vibrational  response  of  cylindrical  shells  to 
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point  excitation.  The  measured  driving  point  admittance  is  compared 
with  predictions  using  the  Mean-Value  Admittance  Method. 

The  scope  of  this  investigation  is  limited  to  circular  cylin¬ 
drical  shells,  although  many  of  the  methods  used,  and  the  results  ob¬ 
tained,  can  be  extended  to  plates  and  shells  of  other  shapes.  The  shells 
that  are  studied  may  be  classified  as  thin;  that  is  to  say  that  the 
stress  distribution  can  be  taken  as  linear  through  the  thickness.  The 
shells  are  homogeneous  and  isotropic.  The  special  case  of  anisotrophy, 
orthotrophy,  is  also  investigated.  Complicating  effects  which  are  not 
addressed  are:  initial  stresses,  variable  thickness,  large  deflections 
(nonlinear  behavior),  and  the  surrounding  acoustic  media.  The  effects 
of  shear  deformation  and  rotary  inertia  are  neglected. 

Quite  often  a  ring  is  added  to  a  cylinder  to  stiffen  the  shell 
at  a  point  of  excitation,  thereby  reducing  the  vibrational  energy  that 
is  transmitted  to  the  shell  structure.  A  radial  point  force  in  the 
plane  of  the  ring  will  excite  one-dimensional  modes  in  the  ring  which 
will  then  couple  with  two-dimensional  modes  in  the  shell.  Analytical 
models  are  evaluated  that  describe  in  a  simplified  manner  the  coupling 
between  shell  and  ring  modes,  while  retaining  the  important  features  of 
the  vibrational  response.  Methods  are  presented  for  the  prediction  of 
the  mean-value  of  the  driving-point  admittance  of  shells,  both  with  and 
without  stiffening  rings.  Although  not  specifically  discussed,  the  re¬ 
sults  of  these  investigations  can  also  be  extended  to  the  use  of  string¬ 
ers  to  stiffen  the  shell  in  the  axial  direction. 

During  the  course  of  this  investigation  numerous  driving- 
point  and  transfer  admittances  were  measured  on  two  cylindrical  shells, 
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one  with  and  one  without  ring-stiffening.  In  order  to  provide  a  complete 
analytical  treatment,  it  was  later  decided  to  limit  the  scope  of  investi¬ 
gations  principally  to  driving-point  admittances.  The  extension  of  the 
mean-value  theory  to  the  calculation  of  transfer  admittances  is  intro¬ 
duced  but  not  pursued. 
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CHAPTER  II 
THEORY 


2.1  Classical  Mode-Sum 

Fundamental  to  most  vibrational  theories  is  the  concept  that 
the  response  of  any  closed  system  to  an  arbirtary  excitation  can  be  rep¬ 
resented  by  the  superposition  of  the  system's  natural  modes: 

5  =  2  5v(x,y,z)  ej(Jt  *  C2.1) 

v 

Here,  Sj  represents  the  scalar  displacement  on  the  surface  of  the  vibra¬ 
tor,  and  are  the  mode  functions.  The  harmonic  time  dependence  of 
eJut  is  assumed  for  all  the  excitation  and  response  functions  and  will 
hereafter  be  omitted.  The  mode  functions  or  eigenfunctions  of  a  system 
represent  natural  states  in  which  the  system  prefers  to  vibrate,  inde¬ 
pendent  of  the  forcing  frequency.  Along  with  satisfying  the  boundary 
conditions  of  the  system,  the  mode  functions  must  obey  the  property  of 
orthogonality.  It  will  suffice  to  say  that  the  property  of  orthogonality 
cannot  be  satisfied  unless  the  system  is  closed;  by  this  we  mean  that  it 
is  reasonably  isolated  from  its  surroundings  such  that  vibrational  energy 
is  not  lost  or  gained  in  an  unpredictable  manner. 

By  representing  the  system  in  terms  of  mode  functions,  the  re¬ 
sponse  of  any  homogeneous  system  can  be  reduced  to  a  set  of  one-degree  of 
freedom  oscillator  equations,  one  equation  for  each  mode.  This  is  accom¬ 
plished  by  incorporating  into  the  elements  of  the  one-degree  of  freedom 
oscillator  information  concerning  the  forcing  function  and  the  point  of 


observation  on  the  surface  of  the  vibrator.  Rather  than  writing  out  the 
oscillator  equations,  it  is  more  convenient  to  use  the  symbolic  notation 
of  electrical  network  theory  to  represent  the  vibrating  system.  The  cir¬ 
cuit  representation  of  a  homogeneous  vibrating  system  is  illustrated  in 
Figure  1.  By  representing  the  vibrator  by  its  electrical  equivalent,  we 
can  make  use  of  the  powerful  theorems  developed  for  electrical  networks 
and  transmission  lines  [11] . 

The  circuit  shown  in  Figure  1  is  called  a  canonic  circuit,  and 
consists  of  the  parallel  connection  of  an  infinite  number  of  series- 
resonant  circuits,  one  for  each  mode.  Each  series  circuit  is  composed 
of  a  mode  mass  (Mj),  compliance  (C^),  and  resistance  (R^,)  .  At  very  low- 
frequencies,  the  system  usually  behaves  like  an  inductance  (M0) ,  which 
represents  the  rigid  body  translation  of  the  vibrator.  At  very  high 
frequencies,  the  system  is  usually  shunted  by  a  capacitance  (C0)  which 
represents  the  compliance  at  the  interface  between  the  driver  and  vibra¬ 
tor,  or  the  local  deformation  of  the  vibrator  material  in  the  vicinity 
of  the  driver. 

Given  the  resonance  frequencies,  ,  and  the  modal  loss  fac¬ 
tors,  nv »  the  mode  compliance  and  resistance  can  be  expressed  in  terms 
of  the  mode  mass: 
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The  general  expression  for  the  mode  mass  is: 
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Here,  we  are  assuming  that  the  surface  mass  density  (m)  is  constant  over 
the  surface  area  (S) .  The  total  mass  (M)  is  then:  M  =  mS.  The  quantity 
<5.j>is  the  average  of  the  square  of  the  yth  mode  function  TOVer  the  sur¬ 
face  area.  The  quantity  is  the  excitation  constant,  and  is  defined 
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The  excitation  constant  states  how  much  of  the  total  force  (FQ)  is  avail¬ 
able  to  excite  a  particular  mode  (Cy).  Term  A  symbolically  denotes  the 
point  of  observation. 


For  a  point  force  the  excitation  constant  becomes: 

Cy(F) 


Cy(A) 


=  1 


(2.6) 


where  F  denotes  the  coordinates  of  the  point  force.  Thus,  the  point 
mode  mass  becomes: 


M  = 


M  <Sy> 

£y(F)£y(A) 


(2.7) 


It  can  easily  be  shown  by  integrating  the  square  of  the  mode  function 

over  the  surface  of  a  vibrator  of  total  mass  M  that  =  — -  for  one- 

di’mensional  vibrators  such  as  beams  and  rings,  and  that  Mv  =  —  for  two- 
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dimensional  vibrators  such  as  plates  and  shells.  Here  we  are  assuming 
that  the  mode  functions  are  sinusoidal,  and  that  the  distortion  of  the 
mode  function  at  the  edge  of  the  vibrator  may  be  neglected. 
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The  admittance  of  the  system  can  now  be  expressed  as  the  sum 


of  the  individual  mode  admittances  (neglecting  MQ  and  CQ) : 
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The  admittance  is  simply  the  resulting  velocity  of  the  system  for  a 
given  force  distribution. 


For  the  driving-point  admittance  (driver  and  receiver  coinci¬ 
dent)  ,  all  of  the  modes  are  excited  with  the  same  phase  as  the  driver, 
and  all  the  circuit  elements  are  positive.  Therefore,  resonances  and 
antiresonances  alternate  with  increasing  frequency  as  predicted  by  Fos¬ 
ters  Theorem.  At  a  resonance  the  admittance  is  a  maximum,  and  at  an 
antiresonance  the  admittance  is  at  a  minimum. 


For  the  transfer  admittance  (driver  and  receiver  separated), 
the  circuit  elements  no  longer  need  be  positive.  The  circuit  elements, 
My,  will  be  negative  if  the  point  of  observation  is  vibrating  in  anti¬ 
phase  with  the  driven  point.  If  the  elements  of  two  successive  reso¬ 
nant  circuits  have  opposite  signs,  their  contributions  add  in  the  fre¬ 
quency  region  between  the  resonance  frequencies.  The  result  is  plotted 
as  a  shallow  trough,  similar  to  the  trough  in  the  transmission  curve  of 
a  band-pass  filter. 


The  mode-sum  procedure  offers  a  straightforward  and  universal 
approach  to  the  description  of  vibrating  systems.  Yet,  it  has  limited 
applications  for  most  of  the  complex  structural  systems  in  use  today. 
The  mode  functions  and  resonance  frequencies  are  usually  too  difficult 
to  predict,  or  if  predictable,  may  be  too  numerous  and  calculations  be¬ 
come  tedious  .  If  only  a  few  resonance  frequencies  are  of  interest,  a 
prudent  approach  for  complex  systems  is  to  experimentally  measure  the 
resonance  frequencies,  either  on  the  actual  system,  or  if  not  feasible, 
on  a  scale  model  of  the  system.  The  procedure  is  further  limited  by  the 
fact  that  there  are  no  closed  systems  in  practice.  The  best  we  can  hope 
for  is  that  the  system  is  only  loosely  coupled  to  its  surroundings,  or 
that  the  exchange  of  vibrational  energy  to  the  system's  surroundings  is 
somewhat  tractable.  Finally,  the  mode-sum  procedure  may  provide  more 
information  than  is  of  interest.  In  fact,  the  detailed  information  pro¬ 
vided  by  the  procedure  may  lead  us  toward  a  false  sense  of  accuracy.  In 
many  practical  structures,  the  resonances  are  so  dense  or  overlapping 
because  of  damping  that  the  identification  of  individual  modal  responses 
would  be  of  little  value.  Instead,  we  are  interested  in  the  gross  re¬ 
sponse  over  some  frequency  band  that  may  contain  many  overlapping  modal 
responses.  A  suitable  measure  of  this  gross  response  would  be  the  geo¬ 
metric  mean  of  the  admittance  over  some  frequency  interval.  This  mea¬ 
sure  is  the  topic  of  the  next  section. 

2. '2  Mean-Value  Oriving-Point  Admittance 

If  the  velocity  amplitude  of  a  vibrating  system  is  plotted  on 
a  logarithmic  scale,  the  resonance  maxima  extend  above  the  mean  line 
through  the  frequency-response  curve  as  much  as  the  minima  extend  below 
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it.  As  the  damping  increases,  the  maxima  and  the  minima  coalesce  to  the 
same  geometric-mean  response  curve.  The  mean  line  can  be  considered  to 
represent  the  characteristic  wave  field  that  would  be  generated  by  the 
driver  if  the  vibrator  were  sufficiently  large  or  sufficiently  damped  so 
that  reflections  from  the  boundaries  do  not  contribute  significantly  to 
the  response  at  the  driving  point. 


By  integrating  the  modal  responses  over  a  unit  frequency  in¬ 
terval,  Skudrzyk  [9,  10,  11]  has  shown  that  the  geometric  mean  of  the 
driving-point  admittance  of  any  homogeneous  vibrator  is  described  by: 


y  =  G  +  jB  =  — - — 
0  0  ° 


(2.9) 


where  e  is  the  frequency  difference  between  successive  resonances  (in¬ 
verse  of  the  density  of  resonances)  and  Mv  is  the  point  mode  mass.  The 
real  part  of  the  admittance  describes  an  energy  carrying  wave  that  prop¬ 
agates  away  from  the  driving  point  (characteristic  wave) ,  and  the  imag¬ 
inary  part  represents  a  wattless  field  that  is  observed  near  the  driver. 


For  vibrators  with  constant  e^M^,  the  imaginary  part  of  the 
admittance  is  zero.  This  result  follows  from  the  fact  that  the  contribu¬ 
tions  to  B0  from  the  modes  whose  resonances  are  above  the  forcing  fre¬ 
quency  cancel  the  contributions  of  the  modes  whose  resonances  are  below 
the  forcing  frequency.  This  applies  to  most  of  the  important  vibrators, 
such  as  rods,  plates,  and  shells  (above  the  radial  resonance  frequency). 
The  exception  is  the  beam  or  ring,  where  evMv  is  proportional  to  uj. 

For  these  types  of  vibrators,  BQ  is  a  negative  mass  admittance  and  is 
equal  to  the  real  part. 
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For  most  vibrators,  the  task  of  determining  the  driving-point 
admittance  essentially  reduces  to  that  of  determining  the  density  of 
resonances . 


CHAPTER  III 


CYLINDRICAL  SHELL 

In  this  chapter,  simple  expressions  are  derived  for  the  pre¬ 
diction  of  the  mean-value  of  the  driving-point  admittance  of  a  circular 
cylindrical  shell .  The  task  of  determing  the  driving-point  admittance  es¬ 
sentially  reduces  to  that  of  determining  the  density  of  resonances.  The 
derivation  of  the  density  of  resonances  is  generalized  to  the  case  of 
orthotropy  where  the  bending  stiffnesses  are  not  equal  along  the  princi¬ 
pal  axes  of  the  material.  The  orthotropic  cylindrical  shell  is  a  useful 
model  of  a  cylinder  with  stiffening  rings  or  stringers.  The  predicted 
mean-value  admittance  of  an  isotropic  shell  is  compared  with  experimental 
measurements,  and  an  interpretation  of  the  data  is  provided.  For  de¬ 
tails  of  the  experimental  investigations,  the  reader  is  referred  to  Ap¬ 
pendix  A. 

3.1  Mode  Parameters 

3.1.1  Mode  Mass 

The  displacement  in  the  radial  direction  of  the  cylindrical 
shell  may  be  described  by  the  mode  functions: 

^m,n  =  ^  cos  n9  cos  km:  »  n-  0,1, 2, 3...  ,  (3.1) 


nnr  a 
=  i  =  a 


wnere: 


k, 


(5.2) 
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The  coordinate  system  and  dimensions  of  the  cylindrical  shell  are  shown 
in  Figure  2.  Here,  n  is  the  integer  number  of  waves  in  the  circumferen¬ 
tial  direction,  and  m  is  the  integer  number  of  half-waves  in  the  axial 
direction.  It  is  convenient  to  introduce  the  normalized  wavenumber, 
o/a,  which  will  be  used  in  the  shell  frequency  equation. 


5.1.2  Density  of  Resonances 

The  density  of  resonances  (or  mode  density)  for  the  cylindri¬ 
cal  shell  has  been  derived  by  Heckl  [7]  and  Skudrzyk  [11].  Heckl  based 
his  analysis  on  a  simplified  equation  for  the  natural  frequencies. 
Skudrzyk,  using  a  less  simplified  form  of  Heckl' s  equation  and  a  unique 
approach,  obtained  more  accurate  results.  Skudrzyk' s  derivation  can, 
however,  be  improved  in  the  region  about  the  radial  resonance  frequency 


of  the  cylinder.  The  following  derivation  uses  Skudrzyk's  approach,  but 
improves  the  accuracy  of  the  results  near  the  radial  resonance  frequency. 
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It  also  generalizes  to  the  case  of  orthotrophy  where  the  bending  stiff¬ 
nesses  are  not  equal  along  the  principal  axis  of  the  material. 

The  resonance  frequencies  of  a  simply  supported  finite  cylin¬ 
der  of  length  i,  radius  a,  and  thickness  h  are  given  by: 


Q2  = 


(l-y2)c4 
(n2  +  a2)2 


joi2  +  a2)2  - 


n2  (4-u) -2-y 
2 (1-U) 


(3.5) 


where  y = h//l2a  is  the  thickness-to-radius  parameter,  y  is  Poisson's 
ratio,  and  (2  =  u)/uq  .  The  frequency  is  normalized  to  the  radial  reso¬ 
nance  frequency  u)Q =  Cp^/a  ,  where  Cp^  is  the  velocity  of  a  longitudinal 
wave  in  a  plate.  Heckl  shows  that  this  equation  is  in  good  agreement 
with  measurements  performed  by  Arnold  and  IVarburton  [2] .  By  using 
y  =  0.3  this  equation  reduces  to: 


n2  = 


0.190 


y2  [(n2  +  a2)2  -2.6  n2  +1.6]  .  (3.6) 


(n2  +  a2)2 

In  order  to  calculate  the  density  of  resonances,  Heckl  approximated  this 
frequency  equation  by: 


n2 


•  T  T~r  +  ^(n2  +a2)1 
L(n2  +c‘)  J 


+  y(n2  +a2) 

while  Skudrzyk  used  the  following  form: 


(3.7) 


Or  * 


(n2  +o2)2 


2^2 


+  y" (n‘  +  o  ) 


(3.3) 


Equations  (3.6)  through  (3.8)  are  plotted  in  Figures  3  through  5,  re¬ 
spectively,  for  Y  3  0.01.  These  figures  display  the  contours  of  con¬ 
stant  frequency  in  wavenumber  space,  and  are  shown  as  continuous  func¬ 
tions  for  convenience  only.  As  shown  in  these  figures,  Equation  (3.S) 
is  a  much  better  approximation  of  the  frequency  equation  than  is 


•  ■■  r,  ■  . 
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Figure  3 


Contours  of  Constant  Frequency  in  Wavenumber  Space  for 
Equation  (3.6)  . 
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Figure  4  Contours  of  Constant  Frequency  in  Wavenumber  Space  for 
Equation  (3.7)  . 


Equation  (3.7),  particularly  for  S1<1  .  However,  neither  approximation 
is  good  in  the  vicinity  of  (2  =  1  .  We  shall  use  Equation  (3.8)  in  our 
calculations,  but  we  must  exercise  caution  in  the  vicinity  of  Q =  1  to 
assure  a  valid  result. 


Recognizing  that  the  bending  stiffness  is  proportional  to 
thickness,  and  thus  y,  we  can  introduce  different  bending  stiffnesses 
along  the  9  and  z  axes  by  introducing  different  y's  for  each  modenumber: 

y.  =  h^/ZTifa  and  y->  =  h-,//l2a  .  (3.9) 


Here,  y ^  is  the  thickness-to-radius  parameter  in  the  circumferential  di¬ 
rection,  and  y^  is  the  thickness-to-radius  parameters  in  the  axial  di¬ 
rection.  Equation  (3.S)  now  becomes: 


=  +  ^V2+  y^a2)2 

(n2  +  cT)2  1  * 


(3.10) 


Equation  (3.10)  can  be  simplified  by  introducing  the  polar  co¬ 


ordinates  : 


:2  =  n2  +  a2 


and  cos*-  4>  =  — j 

r 


(3.11) 


Thus : 


ft2  =  cos2<j>  +  r‘,(y1sin2q  +  y0cos2i)‘ 


=  cos2))  +  r4B2(q) 


(3.12) 


The  relationship  of  these  parameters  to  the  wavenumbers  is  displayed  in 
Figure  6.  Assuming  that  <t>  is  constant,  the  change  in  r  with  respect  to 
frequency  is: 
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3[2  ,  \  aT 

212  —  =  4r  B  (?)  — 
3[2  3(2 


and 


Ar 


A(2  (2  A(2 


A(2 


2r  r‘B2($) 
0=const . 


2r  s(?Vf22-cos\p 


2 rB ( 9 )  cos49 

(22 


•  (»)  *  1  j) 


The  area  between  the  curves  (2  =  constant  and  (2  +  A(2  =  constant  is  given 
by:  tt/2  tt/2 


,  A(2 

A  =  I  r  Ar  dcp  =  ~ 


d£ 


CY-j^sin2^  +  y2cos24>)  /j_  cos4<p 

n2 


(3.14) 


The  lower  limit  of  integration  <j>,  depends  on  the  value  of  n. 
For  n  =  0,  the  frequency  equation  becomes: 

(22  =  1  +  y^o1* 


(3.15) 

The  lowest  frequency  that  is  possible  is  [2=1,  which  corresponds  to 
C  =  0;  this  means  that  the  cylinder  is  infinitely  long.  Therefore,  for 
(2  >>  1  we  must  use  n  =  0  and  =  0.  For  (2  <<  1,  the  lowest  value  of  <p 
occurs  for  n  =  1  (see  Figure  5),  and  the  frequency  equation  becomes: 


ft2  =  7,  T TT~  *  (Y1+Y->tf2)2  •  (3.16) 

(1  +  a  )  1 

Unless  a  is  very  small,  the  membrane  term  a4/(l  +o2)2  =  1  ,  and  2  would 

be  beyond  the  frequency  region  of  interest.  We  can  thus  neglect  Y?c2 
Compared  to  y  .  If  this  were  not  true,  that  is  if  y -,a2  was  the  domi¬ 
nant  term,  Equation  (3.16)  would  reduce  to  the  frequency  equation  of  a 
beam  along  the  c-axis.  Rewriting  Equation  (3.16)  we  have: 
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Qz  =  cos4*?.  +  y* 

and  i 

\  =  cos'1  _  2  =  cos"  V^-aj2)/",2 

a  vE  .  (3.17) 


Here,  cj^  is  the  first  pure  circumferential  resonance  (n  =  1,  c  =  0) . 

Note  that  if  u  is  less  than  ujj  »  sometimes  called  the  cut-off  frequency, 

the  limiting  angle  becomes  imaginary.  The  error  in  the  approximation 

above  is  negligible  as  long  as  a><<u0  but  not  so  low  that  u  -  .  As 

the  frequency  increases  such  that  (1  approaches  1,  the  error  becomes 

larger,  particularly  if  the  shell  is  thick  or  short.  For  (1  =  1,  we  use 

n  =  1  in  an  approximate  form  of  Equation  (3.6): 

Q2  =  1  -  0.91  cos4cpL  +  r4  (y1  sin2t{)L  +  y7  eos2^)  2 

*  0.91  cos4<j).  +  — - -  (Yt  sin2?.  +  y->cos24>i)2 

L  s  in4<>L  1  L  -  L 

Since  is  small,  cos  4>L~  1  and  sin  ^  ;  and  we  have: 

1  2  2 

1  =  o.91  +  -i-(Y +Y2)2  *  0.91  + 

and 

4>l  =  2/y^T  .  (3.  IS) 


The  area  in  the  r,4>  plane  that  corresponds  to  one  resonance  is 
AnAtf*~a/£  ,  and  the  number  of  modes  per  unit  frequency  interval  is 
given  by: 

A 


it/: 


AN  _  AnAc  _  € 

Ail  AQ  2ma 


dt> 


(Y^ sin“j  +  Y7 cos  P)  cos^A 


.  (3.19) 


*/ 

Oi 


v 1-  cos  ^  ~  1  ,  and  the  density  of  resonances  is: 


For  a  »  1, 
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2^ 


(3.20) 


TIius,  at  high  frequencies,  the  mode  density  is  inversely  proportional  to 
the  geometric  mean  of  the  axial  and  circumferential  thicknesses.  If 
h^  =  h-,  =  h,  the  mode  density  reduces  to  /%  LI 2h  ,  which  is  the  same  as 
that  of  a  plate. 


At  frequencies  in  the  vicinity  of  the  radial  resonance  fre¬ 
quency  (ft  =  1),  the  density  of  resonances  becomes: 

rr/2 

AN 

Aft  “  27Ta 

2^ 

+1 


Neglecting  (y^ -yn ) x  compared  to  y^  we  have: 
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AN 


i  1 

2tt  a  2/Ty2 


ln(2/y2) 


/5  ■& 
2v/2  rrh2 


ln(2/Y2) 


=  0.195  —  ln(2/Y2) 


(3.21) 


Thus,  at  frequencies  in  the  vicinity  of  the  radial  resonance  frequency, 
the  mode  density  for  most  practical  purposes  depends  only  on  the  bending 
stiffness  in  the  axial  direction.  This  is  not  surprising  since  a  great 
many  axial  nodes  (n  =  0)  occur  in  the  vicinity  of  oj0  . 


At  very  low  frequencies  (Q<<  1 ) ,  the  angle  <p  approaches  tt/2  and 
we  can  write  <J>  *  y  -  e  .  Substituting  this  <j>  into  Equation  (3.19),  the 
density  of  resonances  becomes: 


This  result  was  obtained  by  Skudrcyk.  By  noting  that  <i>  approaches  7r/2 
as  the  frequency  decreases,  one  is  tempted  to  assume  that  cos4(?  -  0 
and  that  Equation  (3.19)  becomes: 
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it/2 


tt/2 
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AN  L 


d® 
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AH  2TTa  I  (y  sin2<J)  +  Y  cos2p) ,/  cos40 
11  2  ft2 


2rraY1 


-l 

cos 


A! 


-i 

cos 


d<j> 


Ai 


2-na.y 


^  (it/2  -  cos1/?! ) 


l  *  ^  lAi 

-  sin  /ft  =  - 

2iraY,  2  TTaY, 


(3.23) 


This  is  the  underlying  assumption  in  Heckl's  calculations.  However,  the 

ratio  cos--'  cannot  be  considered  zero  within  the  lower  range  of  <p . 
ft2 

Heckl  noted  that  his  experimental  data  was  a  factor  1.1  greater  than  his 
theory.  Thus,  Equation  (3.22)  would  show  excellent  agreement  with 
Heckl's  experimental  data. 

The  density  of  resonances  at  low  frequencies  depends  primarily 
on  the  bending  stiffness  in  the  circumferential  direction.  This  is  in 
contrast  to  the  situation  at  high  frequencies  where  the  density  of  reso¬ 
nances  depends  on  the  geometric  mean  of  the  stiffnesses  in  the  axial  and 
circumferential  direction. 


In  summary,  the  density  of  resonances  (l/em  n)  for  the  ortho¬ 


tropic  cylindrical  shell  may  be  approximated  by: 
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r 

a  l 
40^0, 


w>wo  1 


AN  _  AN  _1_ 
Acj  Aft  (jjq 


a  t 

4  /T  tt  a* 


ln(2/Yj 


_£i_  /ft 

^  4/2  a* 


“<Lo  J 


(3.24) 


Here,  we  have  introduced  the  circumferential  and  axial  stiffness  con¬ 


stants,  defined  by: 


a 


2  _  Vhl 

l  /IT 


and 


Cpf_h2 

/12 


(3.25) 


The  density  of  resonances  for  the  orthotropic  cylindrical  shell  is  il¬ 
lustrated  in  Figure  7,  for  various  values  of  and  cu .  The  curves  in 
this  figure  were  computed  numerically  using  Equation  (3.10)  and  normal¬ 
ized  to  the  high  frequency  mode  density  of  an  isotropic  shell  with  stiff¬ 
ness  constant  a*  .  At  low  frequencies  the  density  of  resonances  for  the 
most  part  is  inversely  proportional  to  the  stiffness  in  the  circumfer¬ 
ential  direction.  If  this  stiffness  is  reduced  by  h  the  density  is 
doubled.  At  high  frequencies  the  density  of  resonances  is  inversely  pro¬ 
portional  to  the  geometric  mean  of  the  circumferential  and  axial  stiff¬ 
nesses.  If  either  one  of  these  stiffnesses  is  reduced  by  ’s  the  density 
is  increased  by  a  factor  of  V2.  If  both  of  these  stiffnesses  are  re¬ 
duced  by  '-i  the  density  of  resonances  is  doubled. 


cnoo 


For  the  isotropic  shell,  a^=a0=as  .  The  density  of  reso¬ 
nances  for  the  isotropic  shell  is  illustrated  in  Figure  8.  This  figure 
is  normalized  to  the  high  frequency  asymptotic  value.  At  low  frequen¬ 
cies  the  mode  density  increases  proportional  to  the  square-root  of  fre¬ 
quency,  just  as  a  membrane  where  extensional  modes  are  predominant.  In 
the  vicinity  of  the  radial  resonance  frequency,  the  mode  density  attains 
a  value  that  varies  inversely  with  the  thickness  to  radius  parameter, y  . 
The  density  peaks  at  a  value  that  is  greater  than  the  high  frequency 
mode  density.  At  high  frequencies  the  shell  behaves  as  a  plate  with  a 
constant  mode  density  of  (aZ/4ctg).  As  shown  in  Figure  3,  the  contours 
of  equal  frequency  approach  that  of  a  quarter  circle  with  its  center  at 
the  origin,  the  same  as  a  plate.  The  density  of  resonances  for  the 
shell  is,  however,  twice  as  great  as  that  of  a  rectangular  plate  of  the 
same  surface  area  and  thickness.  This  is  a  consequence  of  the  degener¬ 
acy  of  the  circumferential  modes  which  need  only  satisfy  the  condition 
of  periodicity  about  the  circumference. 

Strictly  speaking,  the  density  of  resonances  derived  in  this 
section  is  valid  only  for  simply  supported  cylinders.  But,  these  re¬ 
sults  can  be  applied  to  any  cylinder  with  nondissipative  boundaries. 

This  is  because  the  resonance  frequencies  of  a  simply  supported  cylinder 
are  identical  to  those  of  a  cylinder  with  nondissipative  boundaries 
whose  length  has  been  slightly  adjusted  to  match  the  nodal  lines  closest 
to  the  boundary.  As  long  as  the  axial  wavelength  is  small  compared  to 
the  length  of  the  cylinder  [an  assumption  made  in  deriving  Equation 
(3.6)],  the  error  is  negligible. 
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Figure  8  Normalized  Density  of  Resonances  for  an  Isotropic 
Cylindrical  Shell. 
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Mean-Value  Driving-Point  Admittance 


Using  Equations  (2.9)  and  (3.24),  the  real  part  of  the  driving- 
point  admittance  of  an  isotropic  cylindrical  shell,  excited  by  a  radial 
force  at  an  internal  point  is  given  by: 
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where  m3  is  the  mass  per  unit  surface  area  of  the  shell. 

The  density  of  resonances  derived  in  the  previous  section  is 

I?7T 

for  axial  modes  of  the  form  cos  ~z  .  All  of  these  modes  are  fully  ex¬ 
cited  when  the  driver  is  located  on  the  edge  of  the  shell.  However, 
when  the  driver  is  moved  to  an  internal  point  on  the  shell,  not  all 
modes  are  fully  excited.  For  an  internal  driving  point,  it  is  conven¬ 
ient  to  assume  the  driver  is  at  the  center  of  the  shell.  In  this  po¬ 
sition  the  odd  order  modes  (m  =  1,3,5  ...)  are  fully  excited,  and  the 
even  modes  (m  =  2,4,6  ...)  are  not  excited  at  all.  Thus,  the  density  of 
resonances  (1/ev)  at  an  internal  point  on  the  shell  is  half  the  density 
if  the  driver  were  at  the  edge  of  the  shell.  It  follows,  then,  that  the 
driving-point  admittance  at  an  internal  point  on  the  shell  is  half  the 
driving-point  admittance  at  the  edge  of  the  shell.  From  a  physical 
point  of  view,  driving  at  an  internal  point  on  the  shell  means  that  the 
characteristic  wave  is  radiating  into  a  full-plane;  and  we  can  expect 


that  the  velocity  and  admittance  would  be  half  that  of  the  velocity  and 
admittance  if  driving  at  the  edge,  where  the  characteristic  wave  radi¬ 
ates  into  a  half-plane. 

The  imaginary  part  of  the  driving-point  admittance,  Bs,  is 
equal  to  the  real  part,  Gs,  if  u><w0  ,  and  is  approximately  zero  if  u>^0. 
The  magnitude  of  the  driving  point  admittance  is  then: 


For  most  engineering  applications  it  is  sufficient  to  calcu¬ 
late  the  high  frequency  admittance,  20  1og|Ysj  =  -20  log  (8msa|)  ,  and 
the  radial  resonance  frequency,  fQ  =  c^/2ira  .  The  low  frequency  asymp¬ 
tote  of  20  log  j  Ys }  is  then  obtained  by  plotting  a  line  with  a  +3  dB  per 
octave  slope  that  intersects  the  high  frequency  admittance  at  fQ  .  In 
most  cases  the  increase  above  the  high  frequency  admittance  at  fQ  can  be 
neglected. 

Experimental  verification  of  the  foregoing  mean-value  theory 
was  conducted  using  an  aluminum  cylindrical  shell.  The  dimensions  and 
physical  properties  of  the  shell  are  provided  in  Appendix  A.  The  sur¬ 
face  mass  of  the  shell  is  ms  =  p h$  =  21.6  kg/m2  ,  the  stiffness  constant 
is  as  =  Cp£hs//IT  =  12.5  m2/sec  ,  and  the  radial  resonance  frequency  is 
f0  =  Cp£/2:ra  =  5300  Hz  .  Using  Equation  (3.27),  the  high  frequency  mean- 
value  driving-point  admittance  in  decibels  is; 
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20  log  t Y  |  =  -20  log  (8msa?)  =  -66.7  dB  re  sec/kg  . 

The  thickness  to  radius  parameter  is  y  =  hs/a/T2  =  0.009,  and  the  in¬ 
crease  over  the  high  frequency  admittance  at  fQ  is: 

20  log  »  +1.7  dB  re  sec/kg  . 

✓  2  7T 

For  most  engineering  applications  the  increase  of  1.7  dB  can  be  ignored 
because  of  the  narrow  bandwidth  in  which  it  occurs.  However,  in  order 
to  illustrate  some  of  the  obscure  features  of  vibrating  shells,  the 
1.7  dB  increase  at  the  radial  resonance  frequency  will  be  included  in 
our  predictions. 

Figure  9  shows  the  measured  driving-point  admittance  for  both 
undamped  and  damped  cases,  and  the  predicted  mean-value  driving-point 
admittance.  Here  the  driver  is  located  at  an  internal  point,  approxi¬ 
mately  6  cm  from  the  center  of  the  shell.  From  the  undamped  curve  it  may 
be  seen  that  at  low  frequencies  the  density  of  resonances  increases  with 
frequency,  and  at  high  frequencies  the  density  is  relatively  constant. 

In  the  vicinity  of  f0  there  are  a  great  many  axial  modes  that  occur  for 
the  case  where  no  circumferential  modes  exist  (n  =  0) .  The  circumfer¬ 
ential  resonance  frequencies  can  easily  be  calculated  by  letting  a  =  0 
in  the  frequency  Equation  (3.5).  These  resonances  are  marked  by  an  as¬ 
terisk  in  Figure  9.  From  the  undamped  curve  we  see  that  the  resonances 
usually  occur  as  groups  of  axial  resonances  (m  =  1,2,3...)  about  the 
dominant  circumferential  resonances  (n  *  1,2,3...).  Whenever  a  new 
value  of  n  is  possible,  the  axial  resonances  occur  very  close  to  one 
another.  It  is  interesting  to  note  that  as  the  length  of  the  shell  is 


od  Mean-Value  and  Measured  Driving-Point  Admittance  of  Cylindr 
Driven  6  cm  from  the  Axial  Center. 


increased,  the  axial  resonances  would  become  less  pronounced,  and  the 
frequency  curve  would  approach  a  curve  resembling  only  the  circumferen¬ 
tial  resonances.  This  is  because  the  waves  with  larse  axial  components 
must  travel  farther  and  are  attenuated  more  than  the  pure  circumferen¬ 
tial  waves. 

As  damping  is  increased  many  resonances  are  excited  within 
their  bandwidth,  and  the  admittance  approaches  the  predicted  mean-value 
in  Figure  9.  This  is  significant  because  most  cylindrical  structures 
will  exhibit  a  response  similar  to  the  damped  curve  in  Figure  9.  The 
mean-value  admittance  theory  will  then  provide  a  simple  but  accurate 
method  of  predicting  the  mean  response  of  many  practical  cylindrical 
structures . 

Figure  10  shows  the  same  data  as  Figure  9,  except  that  the 
driving  point  is  located  at  the  exact  center  of  the  shell.  In  this  po¬ 
sition  the  driver  can  only  excite  the  odd-ordered  axial  modes.  These 
modes  are,  however,  fully  excited,  and  many  have  resonance  peaks  that 
are  as  strong  as  the  circumferential  resonances.  As  damping  is  increased 
the  admittance  approaches  the  predicted  mean-value. 

Figure  10  illustrates  a  unique  exception  to  the  mean-value 
theory.  Because  of  the  unique  dimensions  of  the  cylinder,  more  of  the 
resonances  in  the  vicinity  of  the  radial  resonance  frequency,  fD,  are 
even-ordered  than  odd-ordered  axial  resonances.  The  result  is  that  3/5 
of  the  resonances  in  the  vicinity  of  the  radial  resonance  frequency  are 
even-ordered  axial  resonances.  These  even  axial  modes,  however,  cannot 
be  excited  by  driving  at  the  center.  Thus ,  instead  of  half  the  mode 


She 1 1  Driven  at  Center. 
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density  for  the  center  driven  cylinder,  we  have  2/5  the  mode  density  in 

the  vicinity  of  fQ.  The  mean-value  characteristic  admittance  at  fQ  is 
2/5 

then  20  log  =  -2  dB  less  than  our  initial  prediction.  If  the  di¬ 

mensions  of  the  shell  are  changed  slightly,  or  if  the  driver  is  moved 
off  center,  the  dip  in  the  vicinity  of  the  radial  resonance  frequency 


would  disappear. 


Figure  11  is  the  same  as  the  damped  curve  in  Figure  10,  except 
that  the  frequency  is  plotted  on  a  logarithmic  rather  than  a  linear 
scale.  The  dependence  of  the  admittance  below  the  radial  resonance 
frequency  is  simply  a  straight  line  with  a  +5  dB  per  octave  slope.  The 
slight  rise  in  the  predicted  mean-value  at  f  is  omitted. 


Predicted  Mean-Value  and  Measured  Driving-Point  Admittance  of  Cylindrical 
Shell  Driven  at  Center  (Plotted  on  a  Logarithmic  Frequency  Scale). 
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CHAPTER  IV 

RING-STIFFENED  CYLINDRICAL  SHELL 


The  vibrational  response  of  a  ring-stiffened  cylindrical  shell 
to  point  excitation  in  the  plane  of  the  ring  is  investigated  in  this 
chapter.  Several  analytical  models  are  developed  to  describe  the  be¬ 
havior  of  the  ring-stiffened  shell.  Each  model  varies  in  its  complexity 
and  its  accuracy  in  predicting  the  driving-point  admittance  in  the  plane 
of  the  ring.  The  model  selection  in  an  engineering  application  would 
depend  on  what  details  of  the  frequency  response  curve  that  the  designer 
is  interested  in.  The  models  presented  in  Sections  4.2.1  and  4.2.2  dif¬ 
fer  primarily  in  the  degree  of  simplification  used  in  describing  the 
coupling  between  ring  and  shell  modes.  Very  simplified  models  can  be 
used  if  one  is  interested  in  either  very  high  or  very  low  frequencies 
relative  to  the  radial  resonance  frequency.  The  high  and  low  frequency 
models  are  presented  in  Section  4.2.3.  For  thick  rings,  the  high  fre¬ 
quency  model  must  be  modified  to  account  for  the  first  ring  thickness 
antiresonance.  This  modification  is  presented  in  Section  4.2.4.  Final¬ 
ly,  the  application  of  ring-stiffening  to  reduce  the  vibrational  vel¬ 
ocity  of  a  cylindrical  shell  at  the  point  of  excitation  is  discussed  in 
.Section  4.3. 

Throughout  this  chapter,  experimental  measurements  are  pre¬ 


sented  to  support  the  various  analytical  models.  For  details  of  the  ex¬ 
perimental  investigations,  the  reader  is  referred  to  Appendix  A. 
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4. 1  Mode  Parameters 

4.1.1  Mode  Mass 

The  displacement  in  the  radial  direction  of  the  ring  may  be 
described  by  the  mode  functions: 


Su  =  £  cos  u0  ,  u  =  1,2,3  ...  ,  (4.1) 

where  u  is  the  integer  number  of  waves  in  the  circumferential  direction. 
The  mode  functions  of  the  cylindrical  shell  are  described  by  Equation 
(5.1).  The  dimensions  of  the  ring-stiffened  cylindrical  shell  are 
shown  in  Figure  12.  The  shell  is  sectioned  axially  in  this  figure  to 
expose  the  details  of  the  ring. 


For  a  radial  point  force  in  the  plane  of  the  ring,  the  general 
expression  for  the  point  mode  mass  of  the  ring-shell  structure  is: 

^ 

Here,  the  excitation  constant,  <v,  is  unity  for  a  point  force.  The  sur¬ 
face  mass  density  is  a  function  of  the  axial  coordinate  as  follows: 


ms  (c) 


M= 


2iraZ 


,  0  <  z  <£. 


Mr 


2ira£  27rrb 


Me 


l^<  z  <  f^+b  \ 


£j+b  <z<l 


(4.3) 


\.  j 

Assuming  that  the  mode  function,  £  ,  is  fairly  constant  across  the  ring 


width,  b: 
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%  =  [m”(£-  b)] 


(Mj.  +  Og  b) 


=  ra”(£  +  b)  +  4  Mr 

=  %  Ms  (1  +X)  +  h  Mr 

=  k  MS  +  h  My 


(4.4) 


In  the  above  expression,  m|'  =  Ms/£.  is  the  mass  per  unit  axial  length  of 
the  shell.  Tne  final  approximation  is  valid  as  long  as  the  width  of  the 
ring  is  much  less  than  the  length  of  the  shell. 


Equation  (4.4)  is  not  surprising,  since  we  would  intuitively 
assume  that  the  point  mode  mass  of  the  total  ring-shell  structure  would 
be  the  sum  of  the  point  mode  masses  of  the  shell  and  ring  alone.  This 
assumption  is  valid,  however,  only  if  the  axial  mode  function  is  rela¬ 
tively  constant  across  the  width  of  the  ring,  and  if  the  length  of  the 
shell  is  much  greater  than  the  width  of  the  ring.  The  first  restriction 
is  equivalent  to  the  statement  that  the  ring  is  a  one-dimensional  vibra¬ 
tor.  The  second  restriction  above  is  equivalent  to  the  statement  that 
the  shell  is  a  two-dimensional  vibrator,  because  if  the  shell  length  is 


reduced  to  the  width  of  the  ring  the  resulting  mode  mass  would  be  the 
same  as  a  one-dimensional  vibrator  with  a  total  mass  of  Mg+  Mr.  From 


:his  argument,  we  see  that  for  a  very  short  cylinder  the  point  mode 

M-  M_ 

iiass  of  the  ring-shell  structure  would  fall  somewhere  between  (  — ) 

4  *2 

M  <  V I 


M,  ♦  Mr 

and  (  - — £  )  . 
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4.1.2  Ring  Density  of  Resonances 

For  bending  vibrations  in  the  plane  of  the  ring,  the  ring  res¬ 
onance  frequencies  are  described  by  [13] : 


•ji  = 

u 


(u  +  1)  [(u  +  1)2-1]  ^r 
/(u  +  1)2  +  1  r2 


1,2,3 


(4.5) 


Here,  r  is  the  radius  of  the  centerline  of  the  ring  and  the  stiffness 
constant  is:  a2  =  Cghr//l2,  where  c^  is  the  velocity  of  sound  in  a  rod. 


The  ring  vibrates  similar  to  a  center-driven  beam  of  the  same 
cross  section  and  length  2rrr,  such  that  the  vibration  is  symmetric  about 
the  driven  point.  At  the  fundamental  resonance  of  such  a  beam,  the  wave 
travels  down  the  beam  and  is  reflected  back  to  the  driving  point  in  phase 
(A  =  2tt r)  .  The  beam  resonance  frequencies  are,  therefore,  given  by: 


<JU  =  2TTU 


cir»4)u  ar  , —  2  ar 

=  2 ITU  -  =  U -  t'uT  =  U2  — 

2mr  r  u  r2 


(4.6) 


This  is  the  asymptotic  value  of  Equation  (4.5).  The  curvature  of  the 
ring  increases  the  effective  stiffness,  unless  the  bending  wavelength  is 
much  smaller  than  the  ring  circumference.  The  first  few  resonance  fre¬ 
quencies  of  the  ring  are,  therefore,  higher  than  the  frequencies  given 
by  Equation  (4.6). 


If  the  ring  is  attached  to  the  shell,  we  would  expect  the  ring 
resonances  to  be  different  than  those  given  by  Equation  (4.5).  The  shell 
will  increase  the  effective  mass  and  stiffness  of  the  ring  at  the  driven 
point.  Because  the  ring  is  much  stiffer  than  the  shell,  we  may  neglect 
the  increase  in  stiffness  created  by  the  shell;  thus  the  primary  effect 
of  the  shell  on  the  ring  is  to  mass  load  the  ring.  The  increase  in 
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effective  mass  at  the  driving  point  is  simply  the  increase  in  the  point 
mode  mass  due  to  the  shell.  This  increase  in  effective  mass  will  de¬ 
crease  the  ring  resonance  frequencies  by  a  factor  of: 


6  = 


Mr  ♦ : 


C4.7) 


Here,  Jjl  is  the  mode  mass  of  the  ring,  and  ^  is  the  mode  mass  of  the 
shell.  The  ring  resonance  frequencies  with  shell  mass  loading  are, 
therefore: 


U) 


u 


(u  +  1)  [(u  +  l)2-  1}  ar 
*^(u  +  1) 2  +  1  r2 


(4.S) 


This  result  is  demonstrated  more  rigorously  in  Section  4.2.1  by  analysis 
of  the  modal  coupling  between  ring  and  shell. 


Experimental  verification  of  Equation  (4.8)  was  conducted  using 
a  steel  ring-stiffened  cylindrical  shell.  Details  of  this  shell  may  be 
found  in  Appendix  A.  Figure  15  compares  the  measured  ring  resonance 
frequencies  with  those  predicted  by  Equation  (4.3).  The  agreement  be¬ 
tween  the  measured  and  predicted  ring  resonance  frequencies  is  good. 


Above  the  first  few  resonances,  the  ring  resonance  frequencies 
may  be  approximated  by: 


(4.9) 


The  change  in  frequency  between  successive  resonances  (inverse  of  the 
density  of  resonances)  is: 
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4^ 

r2  ctj.  /£ 


(4.10) 


The  mode  density  of  the  ring  is,  therefore,  increased  by  a  factor  that 
is  proportional  to  the  fourth-root  of  the  added  mode  mass  of  the  shell. 


4.1.3  Shell  Density  of  Resonances 

The  resonance  frequencies  of  the  orthotropic  cylindrical  shell 
are  given  by  Equation  (3.10).  By  making  the  following  substitutions: 


Yl  “  /l2  r 


CE  hr 


/7I 


Y,  = 


.2  = 


v/n  a 

Cpl  h5 
/12 


a  = 


01 

Ufl 


and 


cp£  _ 


the  shell  frequency  equation  becomes: 


0) 


m,n 


z 


+ 


+ 


(4.11) 


This  is  the  frequency  equation  for  the  orthotropic  shell,  with  the  cir¬ 
cumferential  bending  stiffness  determined  by  the  stiffness  of  the  ring 
and  the  axial  bending  stiffness  determined  by  that  of  the  shell  itself. 
The  membrane  term  remains  essentially  unchanged. 

ar 

Given  that  Equation  (4.11)  reduces  to  =  -j- nz  for  the  pure 

r 

circumferential  resonances  (a  =  0)  and  that  the  ring  resonance  frequen¬ 

cy  2 

cies  are  wu  -  — u  &  ,  one  is  tempted  to  assume  that  the  circumferential 
r‘ 


bending  stiffness  term  in  the  shell  frequency  equation  should  also  be 
modified  by  the  factor  of  g.  This  would  assure  that  the  shell's  circum¬ 
ferential  bending  pattern  would  match  that  of  the  rings  at  the  ring  res¬ 
onance  frequencies.  This  boundary  condition  exists,  however,  only  in 
the  plane  of  the  ring  and  cannot  be  applied  to  the  entire  shell.  The 
mass  loading  factor  g  can  be  applied  to  the  shell's  circumferential 
stiffness  only  if  the  shell  is  stiffened  by  consecutive  rings  closely 
spaced  along  the  entire  shell,  or  if  the  shell’s  response  is  primarily 
dictated  by  circumferential  modes  which  must  match  the  ring  resonance 
frequencies,  such  as  would  occur  near  the  circumferential  mode  cut-off 
frequencies.  These  conditions  are  demonstrated  more  rigorously  in 
Section  4.2.1,  by  analysis  of  the  modal  coupling  between  ring  and  shell. 

It  was  not  possible  to  conclusively  verify  Equation  (4.11)  by 
experimentally  measuring  the  ring-stiffened  shell  resonance  frequencies. 
The  response  of  the  shell  is  composed  of  numerous  and  overlapping  reso¬ 
nance  maxima  which  are  difficult,  if  not  impossible,  to  resolve  (see  Fig- 
H-},. -.Furthermore,  it  is  difficult  to  identify  the  mode  function  as¬ 
sociated  with  a  particular  resonance  maxima.  However,  since  our  goal  is 
to  estimate  the  mode  density,  a  more  reasonable  approach  toward  verifying 
Equation  (4.11)  would  be  to  compare  the  predicted  number  of  shell  reso¬ 
nances  within  a  given  frequency  interval  with  the  measured  number  of 
resonances. 

Table  I  lists  the  shell  resonance  frequencies  predicted  by 
Equation  (4.11).  The  circumferential  resonance  frequencies  (a  =  0)  do 
not  match  the  rings  and,  of  course,  are  not  excited.  The  resonance  fre¬ 
quencies  in  Table  I  are  partitioned  according  to  the  ring  resonance 


Table  I 


Predicted  Resonance  Frequencies  (in  He) 
of  the  Ring-Stiffened  Cylindrical  Shell 
Using  Equation  (4.11). 
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4.2  Mean-Value  Driving-Point  Admittance 

The  ring-stiffened  cylindrical  shell  may  be  represented  by  the 
general  circuit  shown  in  Figure  15.  This  circuit  is  a  reasonable  repre¬ 
sentation  of  the  mechanical  system  as  long  as  the  velocities  at  points 
(1)  and  (2)  are  equal  in  phase  and  magnitude,  or  as  long  as  the  nondis- 
persive  wavelength  in  the  ring  is  much  greater  than  the  thickness  of  the 
ring.  At  very  high  frequencies  the  circuit  representation  must  be  modi¬ 
fied  to  account  for  the  ring  thickness  antiresonance. 

In  Figure  15,  Zr  is  the  impedance  of  the  ring,  and  Zs  is  the 
impedance  of  the  shell.  The  impedances  must  be  modified  to  account  for 
the  modal  coupling  between  the  ring  and  shell.  The  coupling  can  be  de¬ 
scribed  in  various  degrees  of  complexity,  depending  on  what  features  of 
the  vibrational  response  are  considered  important. 

4.2.1  Hybrid  Mode-Sum  and  Mean-Value  Model 

Considerable  physical  insight  can  be  gained  by  the  representa¬ 
tion  of  the  ring  and  shell  by  their  closed  system  mode-sums,  and  by  con¬ 
sidering  how  these  two  component  systems  are  coupled  to  form  a  new 
closed  system.  We  can  expect  that  the- new  system  formed  by  connecting 
the  ring  and  shell  has  resonance  frequencies  that  are  different  than 
those  of  the  ring  and  shell  alone.  The  interaction  between  ring  and 
shell  will  in  general  depend  on  the  degree  of  matching  between  resonance 
frequencies,  on  the  degree  of  similarity  of  their  mode  shapes,  and  on 
how  the  damping  is  distributed  between  the  two  components. 

The  driving-point  admittance  of  the  ring  and  shell  may  be 
represented  as  the  parallel  combination  of  the  ring  modes  in  series  with 
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the  parallel  combination  of  shell  modes  as  shown  in  Figure  16.  Here, 
the  driving-point  admittance  is  VQ  =  v/F,  and  the  mode  parameters  are  de¬ 
fined  by: 
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(4.13) 


uu  are  the  ring  resonances,  <%  are  the  shell's  pure  circumferential 
resonances  (a  =0),  and  are  the  remaining  shell  resonances  (o  t  0) 


The  circuit  representation  in  Figure  16  is  valid  only  if  the 
ring  and  shell  are  attached  at  a  single  point;  that  is,  at  the  driver. 
However,  they  are  attached  along  the  entire  circumference  of  the  ring, 
and  we  must  modify  the  mode  parameters  and  the  circuit  accordingly. 

The  shell  modes  in  Figure  16  are  purposely  arranged  to  illus¬ 
trate  the  coupling  between  the  ring  and  shell.  The  pure  circumferential 
modes  of  the  shell  are  shown  separate  from  the  remaining  shell  modes 
which  have  some  axial  component.  Because  the  shell's  bending  pattern 
must  match  the  rings  bending  pattern  at  the  driving  point,  for  every  ring 
mode  there  is  a  corresponding  shell  circumferntial  mode.  Furthermore, 
since  these  corresponding  pairs  of  ring  and  shell  modes  have  the  same 
bending  wavelength,  then,  n  =  u.  The  remaining  shell  modes  where  o  ?  0 
represent  those  modes  that  will  propagate  vibrational  energy  away  from 
the  ring. 
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With  the  ring  and  shell  attached  we  can  expect  the  new  system 
to  possess  resonance  frequencies  that  are  different  from  the  ring  and 
shell  alone.  Consider  a  single  ring  mode  and  a  single  shell  circumfer¬ 
ential  mode  driven  at  resonance  as  shown  in  Figure  17.  The  resulting 
simple  circuit  illustrates  how  the  ring  and  shell  mode  parameters  com¬ 
bine  to  form  a  new  mode  branch  whose  resonance  frequency  is  different 
from  either  the  ring  or  shell  alone.  Of  particular  interest  is  the  ef¬ 
fective  bending  stiffness  (or  compliance)  of  the  new  system: 


Jeff 


(JjM  +  h  M  ) 
\>  v  r  s' 
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Here,  cov  is  the  new  circumferential  resonance  frequency  that  is  created 
by  coupling  of  the  ring  and  shell.  The  boundary  conditions  require  that 
v  =  u  =  n.  Or  in  other  words,  the  modal  configuration  of  both  the  ring 
and  shell  must  match  the  modal  configuration  of  the  combined  system.  The 
new  circumferential  resonance  frequencies  of  the  combined  system  are 
therefore : 
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(4.15) 


This  is  the  same  result  as  presented  in  Section  4.12.  As  shown  in  this 
section,  the  resonance  frequencies  predicted  using  this  equation  agree 
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well  with  measured  frequencies  on  a  ring- stiffened  shell.  The  approxi¬ 
mation  in  Equation  (4,15)  is  equivalent  to  neglecting  the  stiffness  of 
the  shell  compared  to  that  of  the  ring.  For  the  shell  on  which  the  mea- 
surements  were  performed  the  ring  stiffness  term  (ctf./r  )  is  140  times 
larger  than  the  shell  stiffness  term  (ot4/a4). 

It  must  be  emphasized  that  the  foregoing  analysis  is  only 
valid  for  the  driver  and  receiver  in  the  plane  of  the  ring.  The  shell 
still  possesses  resonances  that  do  not  match  the  ring's  modal  config¬ 
uration  around  the  circumference.  These  shell  resonances  are,  however, 
not  effectively  excited  by  the  ring,  which  acts  indirectly  as  a  driver. 
We  can  detect  the  shell's  own  circumferential  resonances  as  the  re¬ 
ceiver  is  moved  away  from  the  plane  of  the  ring.  This  is  illustrated  by 
Figures  18  and  19.  As  the  receiver  is  moved  away  from  the  ring,  weakly 
excited  shell  resonances  can  be  observed  that  are  not  detected  in  the 
plane  of  the  ring.  These  resonances  are  simply  not  excited  sufficiently 
to  be  observed  at  the  driving-point  on  the  ring.  Of  particular  interest 
is  the  frequency  region  in  the  vicinity  of  the  second  ring  resonance 
(900  Hz).  Here,  we  see  that  much  of  the  vibrational  energy  of  the  sec¬ 
ond  ring  resonance  (Figure  18)  is  coupled  into  other  shell  resonances 
(Figure  19)  as  one  moves  onto  the  interior  surface  of  the  shell. 

This  also  explains  why  the  mode  masses  in  Figure  17  are 

equal  to  ‘JS-  rather  than  for  a  =  0.  The  shell  cannot  be  driven  at  a 
4  2 

purely  one-demensional  vibration  unless  the  ring  stiffness  term  is  equal 

to  the  shell  stiffness  term:  a  /r  =  a  /a.  This  condition  follows  from 

r  s 

Equation  (4.15).  For  most  practical  structures,  this  condition  should 
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never  occur  because  the  ring  is  usually  designed  to  be  much  stiffer 
than  the  shell. 

The  shell  circumferential  modes  which  match  the  bending  pat¬ 
tern  of  the  ring  might  be  more  appropriately  named  "virtual”  modes, 
since  their  effects  are  quantifiable  by  means  of  the  circuit  elements 
Mn,  Cn,  and  Rn,  but  their  existence  is  questionable. 

In  summary,  for  excitation  in  the  plane  of  the  ring  the  pri¬ 
mary  effect  of  the  shell  on  the  ring  is  to  mass  load  the  ring,  and  the 
primary  effect  of  the  ring  on  the  shell  is  to  stiffen  the  shell  in  the 
circumferential  direction.  By  acting  indirectly  as  a  driver,  the  ring 
only  excites  those  shell  modes  that  closely  match  its  bending  wavelength. 
All  other  shell  modes  will  be  weakly  excited.  This  phenomenon  can  be  in¬ 
terpreted  analytically  as  an  effective  increase  in  the  shell's  circum¬ 
ferential  stiffness,  and  the  general  expressions  derived  in  Section 
4.1.5  for  the  orthotropic  shell  apply. 

The  hybrid  mode -sum,  mean-value  representation  of  the  ring- 
stiffened  shell  is  now  presented  in  Figure  20.  The  driving-point  admit¬ 
tance  of  the  new  closed  system  is: 


(4.16) 
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The  admittance  Yra  is  the  mean-value  admittance  of  the  modes  that  corres¬ 
pond  to  a  constant  circumferential  mode  number.  These  axial  type  modes, 
along  with  the  virtual  circumferential  mode  represented  by  Ln,  form  a 
complete  mode  series  for  a  constant  mode  number  v . 


Figure  20  provides  a  good  picture  of  the  principal  loss  mechan¬ 
ism  at  high  frequencies.  The  ^m  can  be  thought  of  as  a  radiation  admit¬ 
tance  which  conducts  vibrational  energy  away  from  the  ring.  A  detailed 
knowledge  of  the  ring  or  shell  loss  factors  is  not  required  as  long  as 
the  shell  is  sufficiently  large  or  damped,  so  that  a  wave  reflected  from 
the  shell  edges  does  not  significantly  contribute  to  the  characteristic 
wave  travelling  outward  from  the  driving  point. 


To  compute  the  value  of  Ym  we  proceed,  as  before,  by  driving 
the  density  of  resonances  em  (for  a  constant  circumferential  mode  number) 
and  substituting  into  the  general  equation  for  the  mean-value  admittance. 


Figure  21  illustrates  the  mode  density  (l/em)  for  the  experi¬ 
mental  ring-stiffened  cylindrical  shell.  This  figure  was  prepared  numer¬ 
ically  by  computing  the  inverse  of  the  frequency  spacing  between  succes¬ 
sive  resonances  and  plotting  this  against  the  geometric  mean  of  the 
resonance  frequencies.  Equation  (4.11)  was  used  for  these  computations. 
A  closed  form  solution  of  these  n  functions  is  so  complicated  that  it 
has  no  practical  use.  A  more  prudent  approach  is  to  approximate  the  de¬ 
sired  mode  densities  within  the  frequency  regions  above  and  below  the 
radial  resonance  frequency. 
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Figure  21  Shell  Density  of  Resonances  for  Constant  Circumferential 
Mode  Number,  n  ,  (Curves  Were  Computed  Numerically). 
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At  high  frequencies  (u)>co0  )  the  shell  frequency  equation  and 
frequency  spacing  between  resonances  can  be  approximated  by: 

ai  a I  a.2 

-m  *  +ir°  -  “n  *  ^ 
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(4.19) 
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Tliis  approximation  is  very  good  for  n  >_  4.  For  small  values  of  n  the 
high  frequency  approximation  can  be  improved  by: 
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This  approximation  is  very  good  for  n  <  4.  At  low  frequencies  (u)<uQ  ) 
the  shell  frequency  equation  can  be  approximated  by: 


which  leads  to: 


“o  0 
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+  cj; 


7  r  V, 

e™  =  — - —  ^  “  »  (4.21) 

nuu0 

where  p2  =  to2 -  w2  .  These  approximations  are  plotted  in  Figure  22  for 
the  frequency  regions  in  which  they  apply.  The  approximate  mode  den¬ 
sities  in  Figure  22  show  good  agreement  with  the  numerically  computed 
mode  densities  in  Figure  21.  Note  that  the  mode  density  for  each  n  ap¬ 
proaches  infinity  at  the  cut-off  frequency,  u)n  =  uv/0  .  Below  the  cut¬ 
off  frequency,  the  mode  density  is  approximately  zero. 


Given  the  foregoing  approximations,  the  mean-value  admittance 
of  the  shell  for  a  constant  circumferential  mode  number,  n  ,  is: 


■x, 


FREQUENCY  (kHz) 


Approximation  of  Shell  Density  of  Resonances  for  Constant 
Circumferential  Mode  Number,  n  . 
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Above  w  the  mode  density  is  inversely  proportional  to  as,  for  ex¬ 
ample,  a  beam.  Thus,  the  imaginary  part  is  a  mass  reactance  and  is 
equal  in  magnitude  to  the  real  part:  Bm  =  -C^.  Below  wo  the  slope  of 
the  mode  density  varies  widely,  and  Bm  is  not  a  simple  function.  Due  to 
the  symmetry  of  the  curves,  Bm  is  approximately  zero  if  averaged  over  the 
entire  low  frequency  region.  Any  approximation  is,  however,  acceptable 
at  low  frequencies  since  the  imaginary  part  of  the  shell  impedance  is 
negligible  compared  to  that  of  the  ring  impedance. 


Below  ojq  we  are  stretching  the  usefulness  of  the  mean-value 
theory.  The  frequency  spacing  between  resonances  is  so  great  that  the 
details  of  individual  resonances  emerge.  In  fact,  the  increase  in  mode 
density  in  the  vicinity  of  the  ring  resonance  frequencies  can  be  inter¬ 
preted  as  a  resonance  rise.  By  equating  the  nt^  mode  impedance  in  Equa¬ 
tion  (4.13)  to  the  mean-value  impedance  [inverse  of  Equation  (4.22)],  we 

can  compute  an  equivalent  ( i)  for  the  single  circumferential  mode.  The 

e 

(i)  for  the  first  three  circumferential  modes  ( a  =  0)  are  shown  super¬ 
imposed  on  the  mode  density  for  a  constant  mode  number,  n  ,  in  Figure  23. 

In  the  circuit  representation  of  Figure  20,  the  circumferential 

modes  are  included  separately,  and  we  must  therefore  terminate  the  axial 

density  of  resonances  at  the  a  =  1  mode  as  shown  in  Figure  23.  Below 

this  first  axial  resonance,  —  is  zero  and  the  impedance  of  the  axial 

em 

branch  approaches  infinity,  thus,  leaving  only  the  impedance  of  the  cir¬ 


cumferential  branch. 


FREQUENCY  (kHz) 


Equivalent  Mode  Density  of  the  First  Three  Circumferential 
Modes  (n  =  l,2  and  3)  Superimposed  on  the  Shell  Density  of 
Resonances  for  a  Constant  Circumferential  Mode  Number,  n  , 
(Mode  Densities  l/em  are  cut-off  at  a =  1) . 


id  Mode-Sum  and  Mean-Value  Prediction,  and  Measured  Driving-Point 
ttance  of  King-Stiffened  Cylindrical  Shell. 
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Figure  25  Hybrid  Mode-Sum  and  Mean-Value  Admittance  Prediction  Using  the  Measured 

Ring  Resonance  Frequencies  and  Loss  Factors  Compared  with  Measured  Drivin, 
Point  Admittance  of  Ring-Stiffened  Cylindrical  Shell. 


We  can  assume  that  most  losses  will  occur  in  the  shell;  thus 
we  let  Hu  =  0  .  The  shell  circumferential  loss  factor,  Hn  ,  is  assumed 
to  be  about  0.02.  These  crude  approximations  are  adequate  for  our  pur¬ 
pose,  since  we  are  concerned  with  the  predicted  mean-value  response 
rather  than  the  resonance  peaks. 


The  mean-value  driving-point  admittance,  predicted  using  the 
model  in  Figure  20,  is  presented  along  with  the  measured  admittance  in 
Figure  24.  As  shown  in  this  figure,  the  hybrid  mode-sum  model  is  an  ex¬ 
cellent  means  of  predicting  the  important  details  of  the  response  curve. 
Figure  25  shows  the  driving-point  admittance  calculated  using  the  mea¬ 
sured  ring  resonance  frequencies  and  loss  factors.  This  figure  clearly 
demonstrates  the  validity  of  the  circuit  representation  of  Figure  20. 


By  modeling  the  ®  t  0  modes  by  their  mean-value  admittance  and 
the  a  =  0  modes  by  their  individual  admittances,  the  mode-sum  computa¬ 
tions  are  greatly  simplified.  There  are  only  eight  ring  resonances  with¬ 
in  0  -  10  kHz ,  and  the  summation  need  only  be  carried  up  to  v  =  10  for 
accurate  results.  In  fact,  because  the  ring  resonances  are  well  sepa¬ 
rated,  excellent  results  can  be  obtained  by  combining  the  contributions 
of  only  the  two  closest  modes. 


4.2.2  Series  Mean-Value  Model 

Returning  to  the  general  circuit  illustrated  in  Figure  15,  we 
see  that  the  impedance  of  the  ring-stiffened  cylindrical  shell,  ZT,  may 
simply  be  represented  as  the  mean  value  of  the  ring  impedance,  Zr,  in 
series  with  the  mean  value  of  the  shell  impedance,  Zs.  The  ring  and 
shell  impedances  must,  however,  be  modified  to  account  for  the  modal 
coupling  between  ring  and  shell. 
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The  primary  effect  of  the  shell  on  the  ring  is  to  mass  load  the 
ring,  thereby  shifting  the  ring  resonance  frequencies  down  by  a  factor 
of  £  that  is  proportional  to  the  square  root  of  the  added. mode  mass  of 
the  shell.  The  spacing  between  resonances  is,  therefore,  modified  by  a 
factor  of  /£,  and  the  ring  impedance  becomes: 
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(4.23) 


where  mr  is  the  mass  per  unit  circumferential  length  of  the  ring. 

The  primary  effect  of  the  ring  on  the  shell  is  to  stiffen  the 
shell  in  the  circumferential  direction.  By  substituting  the  ring's  bend¬ 
ing  stiffness  for  the  shell’s  circumferential  bending  stiffness  in  the 
shell  frequency  equation,  it  has  been  shown  that  the  spacing  between 
resonance  frequencies  may  be  approximated  by  Equation  (4.12).  At  high 
frequencies  (u)>u)q)  the  shell  impedance  then  becomes: 

2(2ev)Mv  _  4crasMs 
^  "  *  "  tt/m  l 

=  8/a/rmsasar  .  (4.24) 

For  most  structures  the  ratio  of  a/r  is  approximately  equal  to  1.  The 
high  frequency  shell  impedance  is  the  same  as  that  of  a  plate  whose  bend¬ 
ing  stiffness  constant  is  equivalent  to  the  geometric  mean  of  the  ring 

and  shell  stiffness  /a, a_  .  At  low  frequencies  (  w < uj  )  the  shell  im- 

^  r  o 

pedance  becomes : 
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(4.25) 


The  low  frequency  shell  impedance  is  primarily  dictated  fay  the  ring- 
stiffened  modes. 


It  should  be  noted  that  Equation  (4.25)  will  become  less  ac¬ 
curate  at  very  low  frequencies,  where  the  shell  response  is  dominated  by 

M 

the  virtual  circumferential  modes,  whose  mode  masses  are  _JL  ,  and  whose 

4  2 

ar 

resonance  frequencies  are  the  same  as  that  of  the  ring:  =  -y  v2  £ • 

These  shell  modes  are  shifted  downward  by  the  factor  of  S  ,  which  would 
make  the  shell  impedance  appear  softer  than  predicted  by  Equation  (4.25). 


Combining  the  complex  impedances  in  the  proper  manner,  produces 
the  expression  for  the  mean-value  driving-point  admittance  in  the  plane 
of  the  ring: 
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T 


(4.26) 


Figure  26  presents  the  predicted  mean  value  and  the  measured 
driving-point  admittance  of  the  ring-stiffened  steel  cylinder.  The  shell 
was  heavily  damped  for  this  measurement.  The  agreement  between  the  pre¬ 
dictions  and  measurements  is  good,  except  for  the  very  low  frequencies. 

The  mean-value  predictions  can  be  improved  at  the  very  low 
frequencies  by  using  the  modified  circuit  representation  illustrated  in 


Series  Mean-Value  Prediction  and  Measured  Driving-Point  Admittance  of 
Ring-Stiffened  Cylindrical  Shell. 
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Figure  27.  To  more  accurately  represent  the  low  frequency  ring-shell 
coupling,  the  shell  impedance  is  represented  by  the  parallel  combination 
of  the  mean  value  of  the  circumferential  mode  impedances,  Zn,  and  the 
mean  value  of  the  noncircumferential  mode  impedances,  Zm. 


Because  the  shell's  bending  pattern  must  match  the  ring's  bend¬ 
ing  pattern  at  the  driving  point,  for  every  ring  mode  there  is  a  virtual 
shell  circumferential  mode.  The  impedance  Zn  is,  therefore,  the  mean 

Me 

value  of  the  mode  impedances  whose  mode  masses  are  and  whose  frequen¬ 
cy  spacing  is  the  same  as  that  of  the  ring: 
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Here,  m^  is  the  mass  per  unit  circumferential  length  of  the  shell. 

The  impedance  of  the  noncircumferential  modes,  Zm,  may  be  ap¬ 
proximated  by  Equations  (4.24)  and  (4.25).  By  using  these  equations,  we 
are  assuming  that  there  are  many  more  noncircumferential  modes  ( a  r  0) 
than  there  are  circumferential  modes  (a  =  0).  Thus,  we  can  remove  the 
a=0  modes  to  form  Zn  without  significantly  changing  Zm.  This  approxi¬ 
mation  is  valid  except  at  very  low  frequencies.  However,  Zm  is  large 
compared  to  Zn  at  very  low  frequencies  and  can  be  neglected.  The  error 
in  using  Equations  (4.24)  and  (4.25)  for  the  impedance  Zm  is,  therefore, 
negligible. 
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Figure  27  Modified  Series  Circuit  Representation  of  the  Ring- 
Stiffened  Cylindrical  Shell  for  Radial  Point  Force 
in  the  Plane  of  the  Ring. 


Figure  23  Modified  Scries  Mean-Value  Prediction  and  Measured  Driving-Point 
Admittance  of  King-Stiffened  Cylindrical  Slid  1 . 
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Combining  the  complex  impedances  in  the  proper  manner  produces 
the  expression  for  the  mean-value  driving-point  admittance  in  the  plane 
of  the  ring: 

Yt  =  - -  .  (4.28) 

1  7  4.  ^m 


Figure  28  presents  the  predicted  mean-value  admittance  com¬ 
puted  using  Equation  (4.28)  and  the  measured  admittance  of  the  ring- 
stiffened  steel  cylinder.  The  agreement  between  theory  and  measurements 
is  good.  The  predicted  mean  value  in  Figure  2S  is  softer  than  the  pre¬ 
dicted  mean  value  in  Figure  26  at  very  low  frequencies,  and  agrees  better 
with  the  measurements. 

4.2.3  High  and  Low  Frequency  Mean-Value  Models 

If  one  is  interested  in  either  the  very  high  or  very  low  fre¬ 
quency  regions  relative  to  the  radial  resonance  frequency,  the  ring- 
stiffened  cylindrical  shell  can  be  represented  by  very  simple  analytical 
models . 


At  low  frequencies  the  response  of  the  ring-stiffened  shell  is 
principally  that  of  the  ring  with  shell  mass  loading  as  shown  in  Figure 
29.  The  resonance  frequencies  are  the  circumferential  resonance  frequen¬ 
cies  of  the  combined  ring-shell  system  as  defined  by  Equation  (4.15)  and 
illustrated  in  Figure  17.  If  the  ring  is  much  stiffer  than  the  shell: 


az 

u)v  =  v2  -I  8 
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(4.29) 


and  the  frequency  spacing  between  resonances  is: 

e  =  & 

V  r 


(4.30) 


The  low  frequency  impedance  is,  therefore: 
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This  is  the  same  as  the  ring  impedance,  except  for  the  factor  of  /6  and 
the  added  mass  loading  of  the  shell,  which  is  approximately  equal  to  half 
the  shell  mass  per  unit  axial  length. 

At  high  frequencies,  the  ring-stiffened  shell  impedance  is 
represented  by  the  series  mean-value  model  as  presented  in  Section  4.2.2. 
The  high  frequency  impedance  is: 

zH  =  Zr  +  zs 

~  mr  ar  /u  (1  +  j)  +  8  ms  as  ar 

=  (8msaras  +  2/g  mr  ar./*T)  +  j  (2/^  mr  ar /w  )  .  (4.32) 

Figure  30  presents  the  mean-value  admittance  using  the  high  and 
low  frequency  models,  and  the  measured  driving-point  admittance  of  the 
ring-stiffened  shell.  Despite  the  simplicity  of  these  models,  they  pro¬ 
vide  a  good  estimate  of  the  mean-value  admittance.  If  the  transition  at 
the  radial  resonance  frequency  is  smoothed,  the  resulting  prediction 
would  show  very  good  agreement  with  the  measurements.  For  most  engi¬ 
neering  applications,  the  high  and  low  frequency  models  are  all  that  are 
required  for  predicting  the  response  of  the  ring-stiffened  shell. 


tnd  Low  Frequency  Mean-Valued  Predictions  and  Measured  Driving 
Admittance  of  King-Stiffened  Cylindrical  Shell. 
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4.2.4  Ring  Thickness  Antiresonance  Model 

For  thick  rings,  the  high  frequency  model  must  be  modified, 
as  shown  in  Figure  31,  to  include  the  effect  of  the  propagation  of  longi¬ 
tudinal  nondispersive  waves  across  the  thickness  of  the  ring.  With  the 
driver  mounted  on  the  inside  of  the  ring,  the  velocity  on  the  outside 
of  the  shell  in-line  with  the  driver  is  not  equal  to  the  velocity  at  the 
driven  point.  This  phenomenon  was  observed  during  the  experimental  in¬ 
vestigations,  and  is  shown  in  Figure  32.  At  very  high  frequencies  (10 
kHz)  the  magnitude  of  the  velocity  on  the  shell  is  approximately  2h  dB 
higher  than  at  the  driving  point.  This  slight  rise  is  due  to  the  first 
resonance  (or  more  correctly,  antiresonance)  across  the  thickness  of  the 
ring. 


The  propatation  of  a  nondispersive  wave  across  the  ring  thick¬ 
ness  can  be  modeled  by  a  transmission  line  (or  rod)  of  length  hr  as 
shown  in  Figure  31  [11,14].  Here,  the  impedances  are: 

khr 

Zj  =  J  Zc  tan  —  (4.33) 

Lm 


and 


sin  kh. 


(4.54) 


where  the  wavenumber,  k  =  w/c^,  and  the  characteristic  impedance  (of  a 
nondispersive  wave),  Zc  =  pc  •  The  ring  material  may  be  considered  in¬ 
finite  in  one  dimension;  thus  the  velocity  of  sound  is: 


V  =  •  '  c4-35) 

As  a  crude  approximation,  the  cross  sectional  area  S  may  be  taken  as  the 
square  of  the  width  of  the  ring,  and  the  characteristic  impedance  is 


then : 


MECHANICAL  SYSTEM 


<Zr 


♦  Zx)  + 


Z2  (zl  ♦  2S) 

Zn  +  Zj  ♦  2S 


CIRCUIT  REPRESENTATION 


Circuit  Representation  of  Ring-Stiffened  Cylindrical  Shell, 
Modified  to  Include  the  Effect  of  Nondispersive  Waves 
Across  the  Ring  Thickness. 


Measured  Velocity  on  Inside  of  King  and  on  Outside  of  Shell  In-Line  with 
Driving  Point,  Driver  is  Mounted  on  Inside  of  King. 
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Zc  =  p  cpl b2  ‘  C4‘36) 

We  are  assuming  that  the  wave  propagation  through  the  ring  thickness  is 
equivalent  to  the  wave  propagation  in  a  rod  of  length  hr  and  cross  sec¬ 
tional  area  of  b2  . 

From  Figure  31  the  ratio  of  the  shell  velocity  over  the  ring 
velocity  (at  the  driving  point)  is: 


This  function  is  plotted  in  Figure  33  for  various  values  of  Zs/Zc.  The 
peak  at  khr  =  ff/2  occurs  when  the  ring  thickness  is  in  antiresonance  at 
the  driving  point.  As  the  shell  impedance  increases  relative  to  the 
characteristic  impedance  of  a  longitudinal  wave  in  the  ring  material, 
the  antiresonance  peak  decreases.  When  Zg  =  ZQ,  the  ring-shell  inter¬ 
face  will  not  reflect  the  incident  wave,  and  vs  =  vr  for  all  frequencies. 


and 


For  the  experimental  ring-stiffened  cylindrical  shell: 

is. .  •"■“A  .  0.3 
Zc  PCp£l>2 

khr  =  0.076  f  (in  kHz)  . 


At  10  kHz,  khr  =  0.76  and  the  shell  velocity  is  about  2\  dB  higher  than 
the  ring  velocity.  This  is  in  good  agreement  with  the  measurements 
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presented  in  Figure  32. 


From  the  equivalent  circuit  in  Figure  31: 


ZT  =  (2r  +  Z1) 


Z2  +  2S  +  Zj 


(4.38) 


Substituting  the  relationships  for  and  ,  the  expression  for  the  high 
frequency  driving-point  admittance  becomes: 


ZH  =  Zr  +  Zs 


(Zc/Zs)  -  j  cot  khj 


(4.39) 


(Zs/Zc)  -  j  cot  khr 

Note  that  if  khT  is  very  small,  this  expression  reduces  to:  Z^  =  Z^  +  Zc 


4.3  Velocity  Insertion-Loss  of  Ring 

Quite  often  a  ring  is  added  to  a  cylindrical  shell  to  stiffen 
the  shell  at  the  point  of  excitation,  thereby  reducing  the  vibrational 
energy  that  is  transmitted  to  the  shell.  A  practical  application  is  the 
mounting  of  a  motor  inside  a  cylinder.  For  a  rotating  imbalance,  the 
motor  may  be  considered  a  constant  force  generator,  and  we  are  interested 
in  reducing  the  velocity  at  the  motor  mountings  by  the  insertion  of  the 
ring.  Defining  the  insertion  loss  as  the  difference  in  decibels  be¬ 
tween  the  shell  velocities  before  and  after  insertion  of  the  ring,  we 
have : 

v2  Z 

IL  =  10  log  \-—\  =  20  log  |— |  dB  .  (4.40) 

V1  Z2 

Here,  the  1  and  2  subscripts  refer  to  before  and  after  insertion  of  the 
ring,  respectively. 


At  low  frequencies  (u<u>0)  the  impedances  above  are: 
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and 


'/ojz 

8ms  as 


=  4/2  ms  a*  /u'J/u  (1  -  j) 


r(l  +  j)  . 


Z.,  =  2/3  (mr  +'inip  arvS 


(4.41) 


(4.42) 


By  combining  these  complex  impedances  in  the  proper  manner  and  taking  the 
magnitude,  the  insertion  loss  is: 


IL  =  20  log 


4  ms  q|  /a>0 
/B  (nij.  +  h  m£)  ar 


-  20  log  u 


(4.43) 


At  high  frequencies  (  oi  >  u)Q )  the  before  and  after  impedances 

are: 

Zj  =  8msa|  (4.44) 

and 

=  8msasar  +  2^^^/!)  (1  +  j)  .  (4.45) 

The  insertion  loss  is,  therefore: 


IL  =  20  log  - - -  =  -10  log[(A  +  B)2  +  A2]  ,  (4.46) 

(A  +  B)  +  jA  . 


where : 


and 


A 


/S  mr 

4  ms 


a. 


a: 


=  a  /a 

r  S 


(4.47) 


(4.48) 


Figure  34  shows  the  mean-value  insertion  loss  computed  for  the 
ring-stiffened  steel  cylinder.  'The  measured  driving-point  admittance  of 
the  damped  shell,  with  and  without  the  stiffening  ring,  is  shown  in 


.  ...» 
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Figure  33.  The  difference  between  these  curves  is  the  velocity  inser¬ 
tion  loss.  The  upper  curve  was  measured  at  the  edge  of  the  shell  and 
was  shifted  downward  6  dB  to  represent  an  interior  point.  The  lower 
mean-value  curve  in  Figure  35  was  computed  by  subtracting  the  insertion 
loss  from  the  upper  mean-value  curve. 

Above  the  radial  resonance  frequency,  the  ring  is  very  effec¬ 
tive  in  reducing  the  velocity  transmitted  to  the  shell.  At  lower  fre¬ 
quencies,  the  shell  is  inherently  stiff  and  the  addition  of  the  ring  is 
not  as  effective.  If  the  shell  was  prestressed,  for  example,  by  ex¬ 
terior  fluid  pressure,  the  shell  would  become  stiffer  and  the  insertion 
loss  of  the  ring  would  be  less  effective. 

By  substituting  mr  =  pbhr  and  a*  =  cEhr/  /12  into  Equation 
(4.46),  one  can  show  that  the  improvement  in  insertion  loss  with  an  in¬ 
crease  in  ring  thickness  and  ring  width  is  approximately: 

AIL  *  -50  1ogAhr  -20  1ogAb  .  (4.49) 

There  is  a  significant  improvement  in  the  insertion  loss  with  an  in¬ 
crease  in  ring  thickness.  However,  with  the  trend  toward  lightweight 
structures,  we  usually  cannot  afford  the  luxury  of  increasing  the  cross- 
sectional  area  or  mass  of  the  ring.  Besides  the  added  mass,  there  is 
another  practical  limit  on  how  much  the  thickness  can  be  increased.  As 
the  ring  thickness  increases,  the  thickness  antiresonance  will  decrease 
toward  the  frequency  region  of  interest,  thereby  increasing  the  shell 
velocity  relative  to  the  driving-point  velocity.  For  example,  if  the 
thickness  of  our  ring  were  increased  by  30  percent.  Equation  (4.49)  gives 
a  lh  dB  increase  in  insertion  loss.  In  Figure  33,  khr  at  10  kilt,  is  now 


Jfc «  >itt 


Measured  Driving-Point  Admittance  of  Cylindrical  Shell  With  and  Without 
Stiffening  Ring. 
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0.76  x  1.3  =  1.0,  and  we  see  that  the  increase  in  shell  velocity  versus 
the  driving-point  velocity  is  3 h  dB.  Tims,  at  10  kHz,  the  benefit  of 
increasing  the  ring  thickness  by  30  percent  would  be  completely  negated 
by  the  effect  of  the  thickness  antiresonance. 

4.4  Transfer  Admittance  in  the  Plane  of  the  Ring 

It  is  not  the  intent  of  this  thesis  to  provide  an  exhaustive 
treatment  of  transfer  admittances.  However,  the  hybrid  model  presented 
in  Section  4.2.1  lends  itself  so  well  to  the  prediction  of  transfer  ad¬ 
mittances  in  the  plane  of  the  ring,  that  it  would  be  wrong  to  pass  over 
this  subject  without  some  discussion. 

We  can  assume  that  the  amplitude  of  the  radial  displacement  in 
the  plane  of  the  ring  is  reasonably  constant  as  the  receiver  is  moved 
around  the  circumference  away  from  the  driving  point  because:  (1)  the 
losses  in  the  ring  are  negligible,  [2)  there  is  no  dispersion  of  the 
vibrational  energy  outward  from  the  driving  point  as  occurs  in  a  two- 
dimensional  vibrator,  and  (3)  there  is  no  edge  distortion  as  occurs  in 
a  finite  beam.  W'ith  receiver  and  driving  point  separate,  but  both  in  the 
plane  of  the  ring,  the  mode  masses  as  shown  in  Figure  20  become: 


Mr  <€u> 
€U(F)  €U(A) 


^r  _1 _ 

2  cos  u <?A 


and 


Ms  <V  > 
b  m.n 


?m,n^  5in,n(^) 


1 _ 

nq>A 


(4.50) 


4  cos 


(4.51) 
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where  $  is  the  location  of  the  receiver,  and  where  v  =  u  =  n.  The  mode 
masses  can  now  be  negative  if  the  point  of  observation  in  in  antiphase 
with  the  driving  point.  Using  Equations  (4.50)  and  (4.51)  in  the  Hybrid 
Mode-Sum  and  Mean-Value  Model,  Figure  20,  we  can  predict  the  mean-value 
transfer  admittance  of  the  ring-stiffened  cylindrical  shell  for  any  lo¬ 
cation  in  the  plane  of  the  ring. 

The  measured  transfer  admittance,  along  with  the  predicted 
mean  value,  is  presented  in  Figures  36  and  37  for  two  separate  receiver 
locations.  The  measured  ring  resonances  and  loss  factors  were  used  in 
these  predictions.  Figure  36  is  with  the  receiver  and  driver  separated 
by  one-half  the  circumference  ($A  =  2tt/2)  .  With  the  receiver  in  this  po¬ 
sition,  cos  v<f)A  =  +1,  -1,  +1,...,  and  the  circumferential  mode  masses 
alternate  in  sign,  producing  shallow  troughs  between  the  ring  resonance 
frequencies.  Figure  37  is  with  the  receiver  and  driver  separated  by  one- 
quarter  the  circumference  (4>A  =  2~/4);  thus,  cos  v4>A  =  +1,  0,  -1,  0,... 
The  odd  ordered  mode  masses  becomes  infinite,  and  the  admittance  is  zero. 
The  response  is,  therefore,  determined  by  only  the  even  ordered  modes, 
whose  mode  masses  alternate  in  sign,  producing  shallow  troughs  between 
their  resonance  frequencies.  This  simply  means  that  for  the  odd-ordered 
modes  the  receiver  is  located  on  a  node  (with  respect  to  the  radial  dis¬ 
placement).  Note,  however,  that  below  the  radial  resonance  frequency 
the  first  and  third  circumferential  modal  responses  are  not  significantly 
suppressed.  The  radial  nodes  associated  with  these  modes  possess  a 
large  tangental  motion.  The  measured  response  is  likely  due  to  the 
transverse  sensitivity  of  the  accelerometer. 


Predicted  and  Measured  Transfer  Admittance  of  Ring-Stiffened  Cylindrical 
Shell  with  Receiver  and  driver  -i  Circumference  Apart  in  the  Plane  of 
the  Ring. 
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Figure  37  Predicted  and  Measured  Transfer  Admittance  of  Ring-Stiffened  Cylindrical 
Shell  with  Receiver  and  Driver  k  Circumference  Apart  in  the  Plane  of 
the  Ring. 
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It  is  significant  to  note  that  the  predicted  transfer  admit- 

ar 

tances  in  Figures  36  and  37  were  computed  using  =  — rn2  B  =  w  in 

2  rz  v 

ar  , 

Equation  (4.11)  rather  than  =  — 5-n2  =  co  /B  >  as  was  used  in  computing 

n  2  v 

the  driving-point  admittance.  The  measured  response  is  much  softer  than 
that  predicted  without  using  the  shift  factor  of  6.  One  may  speculate 
that  for  the  transfer  admittance  in  the  plane  of  the  ring  the  shell's 
response  is  primarily  dictated  by  circumferential  modes  which  must  match 
the  ring  modes.  It  is  not  clear,  however,  why  this  does  not  seem  to  be 
the  case  for  the  driving-point  admittance. 
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CHAPTER  V 

S UNWARY  AND  RECOMMENDATIONS 


S. 1  Summary  of  Results 

Using  the  mean-value  admittance  method  [4],  practical  means 
have  been  demonstrated  for  predicting  the  driving-point  admittance  of 
cylindrical  shells,  both  with  and  without  ring-stiffening.  The  method 
predicts  the  geometric  mean  of  the  response  of  a  vibrator  with  respect 
to  frequency,  and  in  contrast  to  other  methods,  is  applicable  at  low  as 
well  as  high  frequencies.  The  basic  computations  involve  the  mode  mass 
and  the  mode  density.  The  results  obtained  provide  a  gov-d  first  approxi¬ 
mation  with  which  to  evaluate  the  design  of  many  complex  structures  con¬ 
sisting  of  the  fundamental  components  of  plates,  shells,  beams  and  rings. 
Where  the  mode  densities  are  small,  the  mean-value  method  can  be  im¬ 
proved  by  hybrid  mode-sum  representations. 

The  task  of  determining  the  driving-point  admittance  of  the 
cylindrical  shell  essentially  reduces  to  that  of  determining  the  density 
of  resonances.  The  derivation  of  the  density  of  resonances  was  general- 
iced  to  the  case  of  orthotropy,  where  the  bending  stiffnesses  are  not 
equal  along  the  principal  axes  of  the  shell.  The  orthotropic  shell  is  a 
useful  model  of  a  cylinder  with  stiffening  rings  or  stringers.  The  re¬ 
sults  of  the  analytical  development  clearly  demonstrated  that: 

1.  At  frequencies  above  the  radial  resonance  frequency, 
the  effective  bending  stiffness  of  an  orthotropic 
shell  is  equal  to  the  geometric  mean  of  the  bending 
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stiffnesses  in  the  circumferential  and  axial  direc¬ 
tions. 

2.  At  frequencies  below  the  radial  resonance  frequency, 
the  effective  bending  stiffness  of  an  orthotropic 
shell  is  principally  that  of  the  circumferential 
bending  stiffness.  This  result  has  certain  practical 
limitations.  It  is,  of  course,  valid  for  an  isotro¬ 
pic  shell  and  for  a  shell  with  greater  circumferential 
than  axial  stiffness.  However,  if  the  axial  stiffness 
is  much  greater  than  the  circumferential  stiffness, 
this  result  is  valid  at  only  the  very  low  frequencies. 

The  analytical  results  showed  excellent  agreement  with  experimental  mea¬ 
surements  on  an  orthotropic  shell. 

Particular  attention  was  placed  on  the  use  of  ring-stiffening 
to  provide  velocity  insertion  loss  at  the  driving-point.  Several  analyt¬ 
ical  models  were  developed  that  describe  in  a  simplified  manner  the  coup¬ 
ling  between  ring  and  shell  modes,  while  retaining  the  important  features 
of  the  vibrational  response.  Each  model  varies  in  its  complexity  and  its 
accuracy  in  predicting  the  driving -point  admittance  in  the  plane  of  the 
ring.  The  selection  of  which  model  to  use  in  an  application  depends  on 
what  details  of  the  frequency  response  that  the  designer  is  interested 
in. 


The  fundamental  results  of  the  experimental  and  analytical  in¬ 
vestigations  of  the  ring-stiffened  cylindrical  shell  for  point  excitation 
in  the  plane  of  the  ring  are: 


1.  The  primary  effect  of  the  shell  on  the  ring  is  to  mass 
load  the  ring,  thereby  shifting  the  ring  resonance 
frequencies  downward  by  a  factor  that  is  proportional 
to  the  square  root  of  the  added  mode  mass  of  the  shell. 

2.  The  primary  effect  of  the  ring  on  the  shell  is  to 
stiffen  the  shell  in  the  circumferential  direction.  By 
acting  indirectly  as  a  driver,  the  ring  excites  only 
those  shell  modes  that  closely  match  its  bending  wave¬ 
length.  All  other  shell  modes  will  be  weakly  excited. 

This  phenomenon  can  be  interpreted  analytically  as  an 
effective  increase  in  the  shells  circumferential  stiff¬ 
ness,  and  the  general  expressions  derived  for  the  ortho¬ 
tropic  shell  apply. 

The  most  significant  contribution  of  this  investigation  is  a 
simplified  approach  toward  analyzing  the  modal  coupling  between  a  one¬ 
dimensional  and  a  two-dimensional  vibrator  (Section  4.2.1).  This  approach 
uses  the  canonic  circuit  representation  of  a  homogeneous  vibrator.  By 
matching  the  mode  numbers  of  the  two  component  vibrators,  the  mode  pa¬ 
rameters  of  each  vibrator  can  be  combined  by  simple  electrical  network 
theorems  to  produce  effective  mode  parameters  of  the  combined  system. 

The  most  important  analytical  results  of  this  investigation 
are  summarized  in  Table  II. 


5.2  Suggested  Future  Research 


The  principal  shortcoming  of  the  experimental  investigations 


Table  II 


Inverse  of  the  nude  density . 

Sec  section  2.1  for  tfic  general  expression. 


was  that  the  ring-stiffened  shell  was  constructed  with  too  thick  a  ring. 
Consequently,  the  ring  impedance  was  very  large  compared  with  the  shell 
impedance.  Considering  the  accuracy  of  the  driving-point  measurements, 
it  was  difficult  to  observe  the  shell  impedance  which  is  in  series  with 
the  ring  impedance.  More  importantly,  it  would  have  been  desirable  to 


verify  in  more  detail  the  modal  coupling  analysis  presented  in  Section 
4.2.1.  The  fundamental  result  of  this  section  was  the  shift  in  ring 
resonance  frequencies  with  the  addition  of  the  shell: 


Unfortunately,  the  a^/r1*  term  is  140  times  larger  than  the  ct^/a14  term, 
and  the  frequency  shift  factor  3  is  principally  that  of  the  added  mode 
mass  of  the  shell.  In  order  to  verify  Equation  (5.1)  beyond  a  reason¬ 
able  doubt,  it  would  be  desirable  to  test  a  structure  where  the  stiff¬ 
ness  of  the  two-dimensinal  vibrator  could  be  observed  in  the  frequency 
shift.  Such  a  structure  could  most  simply  be  a  rib-stiffened  plate. 
Several  plates  with  varying  rib  stiffeners  could  be  tested  to  verify  the 
shift  in  rib  resonance  frequencies  as  prediced  by  Equation  (5.1)  for 
various  ratios  of  rib  to  plate  stiffnesses.  The  case  where  a(rib)  = 

a(plate)  would  be  of  particular  interest,  because  a  mode  mass  of 

Mnf 

would  probably  have  to  be  substituted  for  -tj—  in  the  frequency  equation. 
Dispite  this  shortcoming,  the  experimental  measurements  did  provide 
strong  evidence  to  support  the  modal  coupling  analysis  in  Section  4.2.1. 


This  is  evident  in  the  Hybrid  Mode-Sum  and  Mean-Value  Model  where  the 
magnitude  of  the  shell  impedance  can  be  observed  between  the  widely 
spaced  ring  resonances. 

As  alluded  to  in  Section  4.3,  a  practical  topic  for  investiga¬ 
tion  would  be  to  study  the  effect  of  various  ring  or  rib  cross  sectional 
geometries  on  the  velocity  insertion  loss.  The  objective  would  be  to  de¬ 
vise  various  methods  of  stiffening  shells  or  plates  to  obtain  an  optimum 
insertion  loss.  Some  examples  of  practical  applications  are  lightweight 
aircraft  and  surface  vehicle  structures,  or  inexpensive  building  par¬ 
titions  . 

The  next  analytical  step  from  this  investigation  is  to  develop 
a  simplified  method  for  the  prediction  of  the  transfer  admittance  of 
cylindrical  shells.  The  transfer  admittance  can  always,  in  concept,  be 
predicted  using  the  mode-sum  procedure.  This  procedure  was  used  in 
Section  4.4  to  compute  the  transfer  admittance  in  the  plane  of  the  ring 
using  the  Hybrid  Mode-Sum  and  Mean-Value  Model.  Here,  there  were  only 
eight  ring  resonances  which  needed  to  be  accounted  for  in  the  calcula¬ 
tions.  Unfortunately,  the  mode  functions  and  resonance  frequencies  are 
usually  too  difficult  to  predict  or  are  too  numerous  such  that  the  mode 
sum  calculation  becomes  tedious.  Furthermore,  for  most  practical  sys¬ 
tems  with  numerous  resonance  frequencies,  we  are  usually  interested  in 
the  gross  response  of  the  system  over  some  frequency  interval,  rather 
than  the  numerous  and  sometimes  misleading  details  of  the  response  pro¬ 
vided  by  a  mode-sum  calculation.  There  is  a  need  for  a  simplified  method 
of  predicting  the  transfer  admittance,  where  we  are  not  interested  in 
the  individual  modal  responses,  but  rather  the  mean-value  of  these 
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responses  over  some  frequency  interval.  Such  a  procedure  was  presented 
by  Williams  [15]  for  vibrating  plates.  Williams  accounts  for  the  spread¬ 
ing  of  the  vibrational  energy  into  two-dimensional  space  and  for  the 
losses  whicli  determine  a  background  or  reverberant  field.  The  results 
presented  by  Williams  for  plates  would,  of  course,  be  valid  for  cylin¬ 
drical  shells  above  the  radial  resonance  frequency.  The  low  frequency 
transfer  admittance,  however,  needs  further  study. 

.Another  analytical  step  from  this  investigation  is  to  consider 
line  or  area  excitation,  such  as  an  acoustical  pressure.  Again,  the 
mode-sum  procedure  can  be  used  by  expanding  the  excitation  in  terms  of 
the  vibrational  modes  of  the  shell,  and  by  computing  a  separate  excita¬ 
tion  constant  and  mode  mass  for  each  shell  mode.  Walter  [20]  has  shown 
that  the  coincidence  of  spinning  acoustic  modes  with  shell  vibrational 
modes  is  a  major, if  not  controlling, factor  in  the  response  of  a  pipe  to 
internal  propagating  sound.  Walter  qualitatively  predicted  the  response 
of  a  pipe  wall  to  a  spinning  acoustic  mode  by  computing  the  excitation 
constant  and  substituting  into  the  mean-value  admittance  expressions. 

His  results,  however,  included  an  effective  length  over  which  the 
acoustic  field  drives  the  shell.  It  is  the  author's  opinion  that  a  re¬ 
lationship  can  be  found  using  the  mean-value  theory  that  is  independent 
if  an  effective  length  of  excitation.  All  of  the  analytical  expressions 
derived  in  this  thesis  are  independent  of  shell  length.  The  mean- 
value  frequency  response  represents  a  characteristic  wave  that  is  inde¬ 
pendent  of  shell  length  or,  in  other  words,  it  represents  a  progressive 
wave  in  the  axial  direction.  Spinning  acoustic  waves  are  also  progress¬ 
ive  in  the  axial  direction.  One  should  be  able  to  evaluate  the  matching 
between  these  waves  on  the  basis  of  a  unit  axial  wavelength  which  for  a 


progressive  wave  is  a  function  of  frequency  and  the  characteristic  vel¬ 
ocity  of  sound  in  the  material.  The  results  of  such  a  study  would  be  of 
great  value  in  designing  pipes  with  high  transmission  loss. 
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APPENDIX  A 

EXPERIMENTAL  INVESTIGATIONS 


A. 1  Test  Specimens 

Various  driving-point  and  transfer  admittances  were  measured 
on  two  circular  cylindrical  shells,  one  with  and  one  without  ring  stiff¬ 
ening.  The  fundamental  dimensions  and  physical  properties  of  these  two 
shells  are  listed  in  Table  III. 


Table  III 


Dimensions  and  Physical  Properties 
of  the  Cylindrical  Shell  Test  Specimens 


Aluminum 

Cylindrical  Sh-11 

Steel  Cylindrical  Shell 
With  Ring-Stiffener 

Shell 

length 

l 

0.724  m 

1.00  m 

nominal  radius 

a 

0.262  m 

0.310  m 

thickness 

hs 

0.008  m 

0.005  m 

Young's  elastic 
modulus 

E 

7.1  x  1010  nt/m2 

1.95  x  1011  nt/mz 

density 

P 

2700  kg/m3 

7700  kg/m3 

Poisson's  ratio 

U 

0.33 

0.28 

...  .. ,  RAng _ 

width 

b 

0.0254  m 

nominal  radius 

r 

0.273  m 

thickness 

h 

0.0655  m 

Young's  elastic 

modulus 

E 

1.95  x  1011  nt/m 

density 

P 

7700  kg/mJ 

Poisson's  ratio 

U 

0.23 
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The  aluminum  cylindrical  shell  was  originally  fabricated  for 
an  underwater  experiment,  and  the  ends  of  the  shell  were  fitted  with 
tongue-and-grooved  rings  for  the  attachment  of  waterproof  caps.  These 
end  rings  present  a  mass-like  loading  on  the  shell  ends,  which  makes  the 
length  of  the  shell  appear  slightly  shorter.  As  long  as  the  driver  is 
located  at  an  internal  point,  the  effective  length  of  the  shell  has 
little  influence  on  the  driving-point  admittance.  The  end  rings,  there¬ 
fore,  will  not  affect  our  measurements. 

Both  the  ring  and  shell  of  the  ring-stiffened  shell  were  fab¬ 
ricated  from  steel,  and  were  attached  by  a  continuous  weld.  The  ring 
was  located  just  off  center  of  the  shell,  so  that  the  ring  would  excite 
all  of  the  shell's  axial  modes. 


A. 2  Method  of  Support 

Both  of  the  cylindrical  shells  were  supported  by  approximately 
six  inches  of  foam  rubber.  This  method  of  support  was  investigated  by 
Williams  in  his  thesis  [15] .  Williams  demonstrated  that  the  vibration  pa~ 
terns  and  natural  frequencies  of  a  plate  were  not  altered  by  the  foam 
support.  The  foam  support  does,  however,  reduce  the  resonance  peak 
heights  if  the  vibrator  is  lightly  damped.  This  is  of  no  importance, 
since  we  are  primarily  concerned  with  the  behavior  of  the  shell  as 
damping  is  increased. 


The  use  of  foam  rubber  is  an  excellent  method  of  supporting 
test  specimens  where  the  objective  is  to  simulate  free-boundary  condi¬ 
tions.  The  impedance  of  the  foam,  although  extremely  small,  is  similar 
to  a  pure  radiation  resistance.  '  The  loss  of  vibrational  energy  to  the 
foam  support  is,  therefore,  predictable  and  constant  for  all  of  the 
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test  specimen's  vibrational  modes. 

A. 3  Method  of  Damping 

Measurements  on  both  shells  were  performed  with  and  without 
damping.  A  commercially  available  compound  was  used  to  dampen  the 
shells.  This  material  is  troweled  on  and  cures  to  a  hard  consistency 
similar  to  plaster.  The  elastic  modulus  and  density  of  the  damping  com¬ 
pound  are: 

E.  =  4  x  109  nt/m2 

and 

pd  =  1600  kg/m3 

In  order  to  compare  the  damped  with  the  undamped  frequency  re¬ 
sponse  curves,  we  must  predict  the  change  in  characteristic  admittance 
of  the  test  specimens  with  the  addition  of  the  damping  compound.  The 
method  used  in  Section  4.2.1  is  used  here  to  estimate  the  influence  of 
the  damping. 


For  convenience,  consider  a  square,  freely  supported  panel  of 
thickness  hp  and  length  t  of  both  sides,  on  which  is  applied  a  layer 
of  damping  compound  of  thickness  hd.  With  the  application  of  the  damping 
layer,  the  mode  parameters  of  the  panel  will  change  as  illustrated  in 
Figure  38.  The  mode  parameters  with  the  1,  2,  and  3  subscripts  are  that 
of  the  panel  alone,  the  damping  layer  alone,  and  the  combined  panel  with 
damping,  respectively.  The  compliance  of  the  combined  system  is: 
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(M£+Mq  . 
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CIRCUIT  REPRESENTATION 


Figure  38  Circuit  Representation  of  Panel  with  Damping  Layer. 
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and  the  resonance  frequencies,  uv3  ,  of  the  combined  panel  with  damping 
is: 


u)2.  = 

Vj 


l" 


(m?  +  n2 


l+nlJMp 


ad 

£4  Cm2 


n\)  Md 


Mp  + 


CA.2) 


Here,  ctp  is  the  plate  stiffness  constant,  and  ad  is  the  stiffness  con¬ 
stant  of  the  damping  layer.  Because  the  modal  shape  of  the  damping  layer 
must  match  the  modal  shape  of  the  panel. 


m->  =  m-i 


n->  =  n-i 


(A.  3) 


and 


uz,  =  —  (m2  +  n?) 

Vj  £4  1  1' 


Mj 


1 


‘V«dMd] 

vMd  J 


=  U) 


vl 


41,  l’l  J 

1  + 

c£  M 

- 2 — E 

.  “p  J 


(A. 4) 


The  change  in  the  resonance  frequencies  of  the  panel  with  the  added  damp¬ 
ing  layer  is,  therefore: 


1  +  5, 


6  * 

u  '!  1  +  6 


where : 


and 


m 


6  ,  fi! 

ap  Mp  Ep  ^p 


<5  »  -  Pd  hd 

m  *  ^  pphp 


J 


(A.5) 


(A.6) 
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The  Ep  and  Ej  are  the  respective  elastic  moduli  of  the  panel  and  damp¬ 
ing  compound;  and  the  pp  and  p^  are  the  respective  densities.  The  6S 
represents  the  effective  change  in  bending  stiffness,  and  5m  represents 
the  effective  change  in  mass. 

Given  the  shift  5^  in  the  resonance  frequencies  of  the  panel, 
the  change  in  frequency  spacing  between  successive  resonances  is  /2C  . 
The  change  in  the  mode  mass  is: 


h  +  h  Mj 


%Mp 


=  1+5 


m 


(A. 7) 


From  the  general  expression  of  the  mean-value  driving-point  admittance. 
Equation  (2.9),  the  change  in  admittance  is  therefore: 


^  Cl  -  6m) 


=  (1  +  5S)  (1+Sm) 


(A.3) 


Tlie  aluminum  cylindrical  shell  was  coated  with  approximately  a 
0.25  cm  thickness  of  damping  compound.  Using  the  shell  data  in  Table  HI. 
the  change  in  effective  stiffness  and  effective  mass  is: 

6S  =  0.0016 


and 

«m  '  0-18 

and  the  change  in  the  characteristic  admittance  is: 

20  1ogAYo  =  -1  dB  re  sec/kg 

This  small  change  was  added  to  the  predictions  in  Section  5.2. 


The  ring-stiffened  cylindrical  shell  was  coated  twice  with  dam¬ 
ping  compound.  The  first  layer  was  approximately  0.15  cm  thick,  and  with 
the  added  layer,  the  final  thickness  was  approximately  0.3  cm  thick.  Tlie 
effect  of  the  first  coat  is: 


6S  =  0.00044 


and 


=  0.058 


20  log  AYq  =  -JjdBre  sec/kg  . 
and  the  effect  of  the  final  coat  is: 


6S  =  0.0039 


and 


5m  =  0.12 


20  1ogAYo  =  -^dBre  sec/kg  . 

Considering  that  the  ring-stiffened  shell  impedance  is  dominated  by  the 
ring  impedance,  these  small  changes  in  the  shell  admittance  were  neglec¬ 
ted. 


Strictly  speaking,  the  estimates  above  are  only  valid  for  fre¬ 
quencies  above  the  radial  resonance  frequencies  where  the  shell  behaves 
as  a  plate.  At  frequencies  below  the  radial  resonance  frequency,  mem¬ 
brane  forces  predominate  and  the  change  in  bending  stiffness  can  be  ne¬ 
glected,  leaving  only  the  change  in  effective  mass.  However,  the  calcu¬ 
lations  above  demonstrate  that  mass  loading  predominates  at  high  fre¬ 
quencies  as  well  as  at  low  frequencies  for  these  particular  shells. 
Therefore,  the  estimates  are  valid  below  as  well  as  above  the  radial 
resonance  frequencies. 

A. 4  Measurement  Apparatus 

The  experimental  measurements  were  condicted  with  the  system 
illustrated  in  Figure  39.  The  measurement  system  was  originally  assem¬ 
bled  by  Hannon  [12]  and  later  improved  by  Williams  [15].  The  basic  prin¬ 
cipal  of  the  system  is  to  maintain  a  constant  driving  force  via  a  feed¬ 
back  loop,  while  harmonically  sweeping  through  the  frequency  spectra. 


L  +  N  2FRA 
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Figure  39  Functional  Diagram  of  the  Measurement  System. 
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For  a  constant  driving  force,  the  measured  velocity  is  directly  propor¬ 
tional  to  the  admittance. 

The  sweep  oscillator  (Spectral  Dynamics  SD104A-S)  generates  a 
harmonic  signal  that  automatically  sweeps  through  a  preset  frequency 
range.  The  output  of  the  oscillator  is  fed  into  an  amplitude  servo/mon¬ 
itor  (Spectral  Dynamics  SD112-1)  which  maintains  a  constant  force  at  the 
driving-point  through  a  feedback  loop.  The  servo  output  is  amplified  by 
a  McIntosh  f  1C- 30  Amplifier  which  drives  the  voice  coil  of  the  Wilcoxon 
Research  F3  Driver.  A  Wilcoxon  Research  Z602  Impedance  Head  containing 
a  force  gauge  and  accelerometer  is  mounted  with  the  driver.  The  force 
and  acceleration  signals  from  the  impedance  head  are  amplified  by  Ithaco 
4S3  Preamps,  and  are  used  by  the  mass  cancellation  unit  to  provide  a 
reference  signal  to  the  servo  that  is  proportional  to  the  force  at  the 
driving-point. 

The  mass  cancellation  unit  was  designed  and  built  by  Williams 
and  the  circuitry  is  described  in  his  thesis.  To  understand  the  purpose 
of  this  component  we  refer  to  part  (b)  of  Figure  40.  Here,  ah  is  the 
acceleration,  is  the  force  measured  by  the  .impedance  head,  is  the 
mass  of  the  mounting  head,  Fs  is  the  actual  force  at  the  driving  point, 
and  Zs  is  the  impedance  of  the  test  specimen.  The  acceleration  at  the 
driving  point  is  the  same  as  the  measured  acceleration;  however,  the 
actual  force  at  the  driving  point  is  related  to  the  measured  force  by: 

F5  ■  Fh  -  "h  “h  ■ 

This  subtraction  is  carried  out  by  an  operational  amplifier  in  the  mass 
cancellation  unit.  The  slight  change  in  the  accelerometer’s  sensitivity 
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Accelerometer 

(measures  a^) 

Driving  Voice 

Coil  (force  generator) 

Force  Gauge 

(measures  F^) 

Mounting  Head  and  Stud 
(combined  mass  of  Mjj) 

Test  Specimen 

(impedance  Zs) 

Figure  40  Impedance  Head:  (a)  Mechanical  System,  (b)  Circuit 
Representation,  and  (c)  Circuit  Representation  with 
Mounting  Compliance. 
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with  frequency  is  also  compensated  for  by  programmable  gains  in  the  oper¬ 
ational  amplifier. 

The  vibrations  of  the  cylinder  are  measured  using  a  1.8  gram 
BBN  Model  501  Accelerometer.  The  accelerometer  signal  is  amplified,  then 
processed  by  a  Spectral  Dynamics  SD101  Dynamic  Analyzer  using  a  10  Hz 
bandwidth  filter  and  RMS  circuitry.  The  output  of  the  RMS  filter  is  con¬ 
verted  from  acceleration  in  decibels  to  velocity  in  decibels  by  subtract¬ 
ing  a  DC  proportional  to  log  frequency  signal  as  follows: 

20  1ogvs  =  20  1ogas  -  20  log  f  -  constant  .  (A.  10) 

The  acceleration  to  velocity  conversion  is  implemented  by  a  simple  volt¬ 
age  divider  network  as  shown  in  Figure  41.  This  circuit  was  designed 
and  built  by  Williams. 

The  dynamic  range  of  the  measurement  system  is  aj^proximately 
60  dB  .  The  adjustment  of  the  system  dynamic  range  involves  a  tradeoff 
between  clipping  the  resonance  peaks  and  losing  the  antiresonances  in 
the  low  frequency  noise  floor.  Since  we  are  seldom  interested  in  the  low- 
frequency  antiresonances,  the  system  was  always  set  up  to  provide  a  true 
plot  of  the  resonance  peaks.  For  some  measurements,  this  means  that  the 
low  frequency  an^iresonances  are  lost  in  the  noise  floor  of  the  system. 
Because  we  are  using  an  accelerometer  to  measure  vibration,  the  noise 
floor  is  proportional  to  acceleration,  and  when  plotted  on  a  velocity 
scale,  will  be  inversely  proportional  to  frequency.  Where  applicable,  the 
svstem  noise  floor  is  indicated  on  the  measurements. 

Two  methods  were  used  to  attach  the  driver  to  the  test  speci¬ 
men.  For  the  ring  stiffened  shell,  the  driver  was  attached  to  the  ring 


i'r  VmT  litr’T'T  -  '* 
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Figure  41  Electrical  Circuit  for  the  Conversion  of  the  Acceleration 
Signal  to  a  Velocity  Signal. 
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by  a  stud  which  was  screwed  into  the  impedance  head  and  ring.  For 
all  of  the  other  measurements,  the  driver  was  attached  to  the  test  speci¬ 
men  by  a  flat  mounting  stud  which  was  glued  to  the  shell  and  screwed  into 
the  driver. 


.Any  method  of  attachment  will  produce  a  small  compliance  be¬ 
tween  the  impedance  head  and  the  test  specimen  [16].  Even  if  the  impe¬ 
dance  head  were  ideally  attached  to  the  test  specimen,  local  deformation 
of  the  test  specimen  material  or  deformation  of  the  impedance  head  itself 

would  produce  a  compliant  effect  at  high  frequencies  [17].  For  the  Z602 
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impedance  head  this  compliance  is  on  the  order  of  1 x 10  m/nt.  This 
compliance  can  usually  be  neglected  compared  to  the  compliance  of  a  glue 
layer  or  other  method  of  attachment.  The  effect  of  a  compliance  between 
the  impedance  head  and  test  specimen  is  illustrated  by  the  equivalent 
circuit  in  part  (c)  of  Figure  40.  In  this  figure,  C  is  the  compliance 
due  to  the  attachment  method,  and  vg  is  the  actual  velocity  of  the  test 
specimen  at  the  point  of  attachment.  From  Figure  40: 

zs  J~ 

Fs  =  zsVs  •  — ^  vh 
Zs+  ju>C 


and 


1  +  juC  Z$ 


20  log  |^-|  =  20  log  |l  +  jo»CZs| 


(A.ll) 


At  high  frequencies  the  compliance  short  circuits  the  driving  velocity, 
thereby  reducing  the  velocity  of  the  test  specimen.  If  Zs  is  large,  the 
measured  acceleration  (a^)  will  be  greater  than  the  actual  acceleration 
(as)  of  the  test  specimen. 

The  compliance  of  both  methods  of  attachment  were  measured  by 

Williams.  The  combined  compliance  of  the  glue  layer  and  flat  mounting 
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stud  was  found  to  be  approximately  1 x 10  m/nt,  and  the  compliance  of 

the  stud  screwed  into  both  the  head  and  test  specimen  was  approximately 
-9 

7x10  m/nt.  The  largest  impedance  measured  with  the  glued  mounting 

4 

was  approximately  95  dB  (  Zs  =  5.6x10  kg/sec  ).  Thus,  the  difference 
between  a^  and  ag  at  10  kHz  is  approximately: 

a. 

20  log  =  31  dB 
s 

The  largest  impedance  measured  with  the  stud  screwed  into  both  the  head 
and  test  specimen  was  approximately  120  dB  (  Zg  =  1  x  106  kg/sec  ) ,  and 
the  difference  between  a^  and  as  at  10  kHz  is  approximately  53  dB.  In 
summary,  at  high  frequencies  and  for  large  test  specimen  impedances,  the 
acceleration  measured  by  the  impedance  head  (a^)  is  grossly  in  error  of 
the  true  acceleration  of  the  test  specimen  (ag) . 

The  problem  of  the  mounting  compliance  was  avoided  by  using  a 
separate  accelerometer  to  measure  the  response  of  the  cylindrical  shells. 
The  impedance  head  acceleration  signal  was,  however,  used  by  the  mass 
cancellation  unit  because  it  correctly  represents  the  acceleration  of  the 
head  mass  (M^) . 

The  attachment  of  the  accelerometer  to  the  test  specimen  can 
influence  the  measurement  in  two  ways.  First,  if  the  mass  impedance  of 


the  accelerometer  approaches  the  impedance  of  the  test  specimen,  the 
added  mass  of  the  accelerometer  will  reduce  the  response  of  the  test 
specimen  at  the  point  of  attachment.  Problems  with  transducer  mass  load¬ 
ing  are  well  known  and  are  discussed  throughout  the  literature  [18,19]. 
Second,  any  method  of  attachment  will  result  in  a  small  compliance  be¬ 
tween  the  accelerometer  and  test  specimen,  and  the  mass  of  the  acceler¬ 
ometer  will  vibrate  in  antiresonance  with  the  compliance  at  the  point  of 
attachment.  Both  of  these  effects  are  illustrated  by  the  equivalent 
circuit  in  Figure  42. 


From  Figure  42,  we  see  that: 


Zs  +  Zl 

vo  =  — -  V1 

Ls 


(A. 12) 


(A. 13) 


where : 


juM  +  -t— r 
jwC 


1  -  w2M  C 


(A. 14) 


The  difference  in  decibels  between  the  measured  velocity  (v->)  and  the 
actual  velocity  (vQ)  of  the  test  specimen  without  the  accelerometer  mass 


loading  is,  therefore: 


20  log  1^— I  =  -20  log 

vo  to2  4s 


(A. 15) 


The  u2/u2  term  represents  the  effect  of  the  antiresonance  between  the 
mass  of  the  accelerometer  and  the  compliance  at  the  point  of  attachment. 
The  antiresonance  frequency  is  defined  as  i j2  =  1/MC  .  The  jwM/Zs  term 
represents  the  effect  of  the  accelerometer  mass  loading  on  the  test 


specimen.  The  effect  of  the  mass  loading  reduces  the  measured  velocity, 
while  the  effect  of  the  antiresonance  between  M  and  C  increases  the  mea¬ 
sured  velocity.  Because  of  their  frequency  dependence,  both  these  ef¬ 
fects  are  only  important  at  high  frequencies. 

The  accelerometer  was  attached  to  the  cylindrical  shells  by 

means  of  a  sticky  putty,  known  commercially  as  Duxseal.  The  compliance 

of  a  small  layer  of  putty  was  measured  by  Williams,  and  was  found  to  be 
-8  , 

approximately  2x10  m/nt.  The  mass  of  the  accelerometer  is  1.8  grams. 
The  effect  of  the  antiresonance  between  M  and  C  at  the  highest  frequency 
of  interest,  10  kHz,  is  then  =  0.14  ,  which  may  be  neglected  com¬ 

pared  to  1  in  Equation  (A. 15).  At  10  kHz,  the  magnitude  of  ju)M/Zs 
should  be  less  than  0.25  for  its  contribution  to  be  neglected  in  Equation 
(A. IS).  This  means  that  Zs  should  be  greater  than  452  kg/sec  (55  dB)  or 
that  Ys  should  be  less  than  -55  dB  .  The  measured  resonance  peaks  of  the 
undamped,  unstiffened  shell  exceed  an  admittance  of  -55  dB  and  are 
slightly  compressed  by  the  accelerometer  mass  loading.  Fortunately,  we  are 
not  as  concerned  about  the  resonance  peaks  of  the  undamped  measurements 
as  those  of  the  damped  measurements.  The  damped,  unstiffened  shell  and 
all  of  the  ring  stiffened  shell  measurements  do  not  exceed  the  limiting 
value  of  -55  dB  and  are  not  affected  by  the  accelerometer  mass  loading. 
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